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Two types of parameter dependent generalizations of classical matrix ensembles are defined 
by their probability density functions (PDFs). As the parameter is varied, one interpolates 
between the eigenvalue PDF for the superposition of two classical ensembles with orthogonal 
symmetry and the eigenvalue PDF for a single classical ensemble with unitary symmetry, 
while the other interpolates between a classical ensemble with orthogonal symmetry and 
a classical ensemble with symplectic symmetry. We give interpretations of these PDFs in 
terms of probabilities associated to the continuous Robinson-Schensted-Knuth correspon¬ 
dence between matrices, with entries chosen from certain exponential distributions, and 
non-intersecting lattice paths, and in the course of this probability measures on partitions 
and pairs of partitions are identified. The latter are generalized by using Macdonald polyno¬ 
mial theory, and a particular continuum limit — the Jacobi limit — of the resulting measures 
is shown to give PDFs related to those appearing in the work of Anderson on the Selberg 
integral. By interpreting Anderson’s work as giving the PDF for the zeros of a certain ratio¬ 
nal function, it is then possible to identify random matrices whose eigenvalue PDFs realize 
the original parameter dependent PDFs. This line of theory allows sampling of the original 
parameter dependent PDFs, their Anderson-type generalizations and associated marginal 
distributions, from the zeros of certain polynomials defined in terms of random three term 
recurrences. 


1 Introduction 

This paper is a companion to our work |p^ . In |p^ we studied the correlation functions for the probability 
density functions (PDFs) 
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where C is the normalization (throughout C will be used to denote some normalization; we remark too 
that it is required A < 1 for (O) to be normalizable) and 

Xi > X2 > ■ ■ ■ > X2n > 0. (1.3) 

The first of these was isolated [Q in the context of a study of generalizations of Ulam’s problem — 
the computation of the distribution of the length of the longest increasing subsequence in a random 
permutation. The second, although not given explicitly in the same paper that ( 0 ) was noted, has its 
origin in a particular model introduced in Q. 
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In we also studied the correlation functions for the particular parameter dependent PDFs 
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where the condition A < a + 1 is required for (]^) and (0) to be normalizable and 

1 > Xi > a;2 > • • • > X2n > 0. 

(In 1^ these PDFs were defined on (—1,1) rather than (0,1) as done here; to define the former simply 
change variables Xj i—> (1 — X2n+i-j)/‘2. in the above.) The PDFs 0 and ( |l.5[ ) were identified in 
[p^ as the only parameter dependent extensions of classical matrix ensembles with an even number of 
eigenvalues, in addition to 0 and 0, with the special property that after integrating over every 
second eigenvalue the eigenvalue PDF of a matrix ensemble with symplectic and unitary symmetry 
respectively results. As noted in by scaling the variables and parameters 
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and taking the limit L —>■ oo, (1.4) and ( |1.5| ) reduce to ( |l.lD and ( |l.2D respectively. 

Let us remark at this point how the above PDFs relate to the classical matrix ensembles. Following 
the notation of [Q, we specify a matrix ensemble with orthogonal {(3 = 1), unitary (/? = 2) or symplectic 
{f3 = 4) symmetry by the eigenvalue PDFs 
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A classical matrix ensemble then refers to an eigenvalue PDF of the form (1.7) with g{x) a classical 


weight function — Gaussian (e ), Laguerre (x"e Jacobi (x"(l —x)^) or Cauchy ((1 + x^) "). Thus 
we recognize 0) in the case A = 0 as the Laguerre orthogonal ensemble with a = 0 and 2n eigenvalues 


which we denote as LOE 2 n|a= 0 ; while we recognize (1.4) in the case A = 0 as the Jacobi orthogonal 
ensemble with a (a — l)/2, 6 = 0 and 2n eigenvalues which we denote as JOE 2 ra| a=(a-i]/ 2 . In the limit 


A -s- —oo of (1.1), (1.4) each pair of coordinates collapses, and after appropriate rescaling and renaming 
of the collapsed pairs the limiting PDFs are 
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The first of these is the Laguerre symplectic ensemble with parameter a = 0, denoted as LSE„|a=o, 
while the second is the Jacobi symplectic ensemble with a a + 1, 6 = 0, denoted as JSE„|a«a+i. 
Consequently we have that ( 0 ) and © interpolate between particular orthogonal ensembles and 
symplectic ensembles. Regarding the PDFs (|l.2| ) and (|l.5D , we require the fact that superimposing 
two orthogonal ensembles as specified by (|1.7D with /3 = 1 at random gives the eigenvalue PDF 
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and ( |l.5|) are recognized as the superimposed ensembles LOE„|a=o U LOE„|a=o 

—oo limit (O) and (|1.5D effectively reduce 
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respectively, which specify the Laguerre unitary ensemble with a = 0, denoted LUE„|a=o, and the Jacobi 
unitary ensemble with b = 0, denoted JUEn,|b=o- Consequently (ED and (ED interpolate between 
particular superimposed orthogonal ensembles and unitary ensembles. 

Our objective in this paper is to give interpretations of each of 0, (y), 0 and (ED as proba¬ 
bility densities relating to longest increasing subsequence problems in settings analogous to that already 


known for (1.1), and also to specify parameter dependent random matrices which have these distributions 
as their eigenvalue PDF. We will see that these two pursuits are intimately related via a conditional PDF 
to be referred to as the Anderson density. To establish this link requires first generalizing the most imme¬ 
diate interpretation of the PDFs in the setting of longest increasing subsequence problems/ last passage 
percolation, discussed in Section 2, to a measure on pairs of partitions A, k with X/k a horizontal strip 
suggested by Macdonald polynomial theory. This is done in Section 3. Taking a particular continuum 
limit, already known from Section 2 as the Jacobi limit, gives rise to the Anderson density. The crucial 
point in making the link with random matrix theory is an interpretation of the workings of Anderson’s 
paper Q as giving the density of zeros of a certain random rational function. This same random rational 
function occurs as part of the characteristic equation for the random projection of a fixed matrix with in 
general degenerate eigenvalues. Such a random projection is used in Section 4.1 to give the construction 
of random matrices with eigenvalue PDFs (El and (El)- It is shown in Section 4.2 that the intricacies 
of the relationship between the Anderson density and the random rational function allow a generalization 
of the joint densities ( |1.4D and (ED — as well as an associated marginal density — to be sampled from 
the zeros of a polynomial generated by a random three term recurrence. 

The Anderson density and the corresponding random rational function have a well defined Laguerre 
limit giving rise to a generalization of the joint densities (ED and ( |l.2| ). As noted in Section 5.1 the 
Laguerre limit of the random rational function occurs as part of the eigenvalue equation for a random 
rank 1 projection of a fixed matrix with in general degenerate eigenvalues, allowing for the construction 
of random matrices with eigenvalue PDFs (ED and (ED- Furthermore it is shown in Section 5.2 that the 
Laguerre limit of the random three term recurrences of Section 4.2 generates polynomials from which a 
generalization of (Eff| ) and (ED , and an associated marginal density corresponding to the Laguerre limit 
of the Selberg integral, can be sampled. 

In Section 6 we carry out a further limiting analysis of the results of Section 4, this time characterized 
by the Jacobi type weights associated to the Anderson density degenerating to Gaussian weights. A 
special case of the resulting parameter dependent Gaussian ensemble is known from Section 2.3 as the 
joint probability density for a particular last passage percolation model. 

In the course of our study we encounter the need to further develop/ reformulate some existing theory, 
in particular the Robinson-Schensted-Knuth correspondence. We also encounter some consequences of our 
findings by way of new insights into some existing results. Such points are presented in the Appendices. 
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2 Last passage percolation and tableaux coordinates 


The generalizations of Ulam’s problem of interest to us can in turn be regarded as generalizations of 
a last passage percolation model introduced by Johansson To dehne the latter consider the right 
quadrant square lattice {{i,j) ■ i,j G Z+}. Associate with each lattice site (i,j) a random non-negative 
integer variable Xij^ chosen from the geometric distribution with parameter aibj so that 

Pr{xij = fc) = (1 - aibj){aibj)^. (2.1) 


For given non-negative parameters oi, 02 , •. •, 6 i, 62 , ■ • ■ the quantity of interest is the distribution of the 
so-called last passage time 

L(ni, 712 ) := max Xi^j ( 2 - 2 ) 

(l,l)u/rh(r!,i ,rt 2 ) 

where the notation (1, l)u/rh(ni, 712 ) denotes that the sum is over all lattice points in a path starting 
at (1,1) and finishing at ( 711 , 712 ), with segments which are either up or right horizontal (such a path 
is said to be weakly increasing). The celebrated Robinson-Schensted-Knuth (RSK) correspondence (see 
e.g. p0| , 1^ ) gives a bijection between tii x 772 non-negative integer matrices and pairs of semi-standard 
tableaux of the same shape /r = (/ii, /72,..., /7„) say with the crucial feature that = L(77i, 772 ). It is 
these variables which after rescaling give rise to the parameter dependent PDFs listed in the Introduction. 

We require some (mostly) known facts about the RSK correspondence in the situation that there is 
a measure on the space of non-negative integer matrices as implied by Q- First, the probability an 
77i X 772 non-negative integer matrix with this measure corresponds to a pair of semi-standard tableaux 
with shape one of content 771 , the other of content 772 is given by 

ni 712 

nn (1 - a^bj)Sf^{ai )Sm(&1 ,..., bji^ ), (2.3) 

i=li=l 

where denotes the Schur polynomial. Second, for 772 > 77i, the joint probability that an 771 x (772 -I-1) 
non-negative integer matrix with measure implied by ( ^.l| ) corresponds to a pair of semi-standard tableaux 
with shape n, content 771 and 772 -f 1, and that the 771 x 772 bottom left sub-block corresponds to a pair 
of semi-standard tableaux with shape k, content 771 and 772 is 


ni n2 + l 

n Y{i^-0'2bj)s^{ai,...,ani)si^{bi,...,bn2)bnJ+i^' 

i=l j=l 


(2.4) 


where 


^1 > Ki > ^2 > K2 > • • • > Mm > > 0. 


(2.5) 


Note that by the assumption 771 < 772 , we require £{^) — the number of non-zero parts of m — to be less 
than or equal to 77i for (^^) to be non-zero. If instead 771 > 772 , then (^) holds with 




( 2 . 6 ) 


and (2.5) must be modihed to read 


Ml > Ki > /i2 > K2 > • • • > Mm > > Mm-l-i > 0- 


(2.7) 


Third, in the case of matrices symmetric about i = j (here i denotes the row counted from the bottom), 
when the RSK correspondence maps the matrix to a single semi-standard tableau, with the parameters 
of the geometric distribution chosen so that 

Pt{x^j = k) = {1 - a^aj){a^aj)^, i < j Pr(a;j_i = fc) = (1 - ai)af, (2.8) 
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the derivation of (2.3) given in 1^ implies the probability that the tableau has shape fi is given by 
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Fourth, in this latter situation, the correspondence is such that 








( 2 . 10 ) 


i=i 


i=i 


In 


(^.10) is given in the form 


i(^,) 


( 2 . 11 ) 

1=1 1=1 

where denotes the length of the jth column of the diagram of /i, which is the sum of the number of 


odd columns; simple reasoning shows that (2.11) is equivalent to (2.10). In Appendix A we will give a 
self contained derivation of (2.10|) which is in keeping with the interpretation of the RSK correspondence 


as a cascade of growth models given in ||^. The same ideas will be used to derive (2.4). 

It follows from (|2.9|) and (2.10|) that if (^.8[) is modified so that the second equation reads 


Pr(a;i,i = fc) = (1 - aai){aai)^ 

then the probability that the tableau has shape /r is given by 

n 

]^(l-aaj) (I - s^(ai,..., o„). 

2=1 


( 2 . 12 ) 


(2.13) 


We will show that (1.1) and (|1.4D result from this probability. The PDFs (1.2) and (|1.5|) will be derived 


from (2.4). 


Before undertaking the derivations, we note that from the origin of (2.4) as a joint probability with 
the corresponding marginal density given or implied by (|2.3| ), the former must satisfy special identities 
with respect to summation over k and summation over /i. Thus let R denote the region (|2.5|) in the case 


712 > n-i, and the region (^) in the case ni > 712 . Then since summing (^) over k G R must give (2.3) 
with 772 1 -^ 772 + 1, it follows 


E ■s«(6 i,..., = s^(&i,..., 6„,+i) 


(2.14) 


k:k^R 


where \fJ-\ = similarly the meaning of |«;|. In fact (2.14) is a well known recurrence satisfied 


by the Schur polynomials special case of (5.10) pg. 72]. Similarly, summing (^) over fi must give 
(p.3[) with K and so 


J^(l (7j5yj2 + l ) ^ ■ ■ ■ : ^ni }^n2 + l ■ i ^ni ) 


(2.15) 


i=l 


which is also a known Schur polynomial identity |^, special case of (1) pg. 93]. 


2.1 Jacobi limit 


The parameter dependent Jacobi ensembles (O), (El are obtained from ( |2.13D , ( |2.4| ) respectively by 
specializing the parameters {ui}, {bj}. In (|2.13D we choose 


{ai,... ,an) = {z, zt,zf,... ,zt'^ ^) 


(2.16) 
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while in (2.4) we choose 

(ai,...,a„J = = {z, zt, zt^,..., zr‘^~'^). (2.17) 

The probabilities then assume an explicit form in terms of the parts of p. and k due to the evaluation 
formula iS 


SA 


___ 1 4-\i — \j—i+j 

(i.(.c-‘) = lEi.c-w. n 

l<2<j;<n 




l<2<j<n 


(2.18) 


where in the second line, which follows from the first by simple manipulation, := (1 — t)(l — 

^2)... (1 _ j jg arbitrary. 

Substituting (2.18) with n* = n in (2.13|) we deduce the following result. 


Proposition 1. On each site {i,j) ofthenxn square lattice, specify a non-negative integer Xi^j according 
to the probability distribution 

Pr(x*,j =k) = (1 - z2f+J-2)(22^z+j-2)fe . ^ J 

= k) = {1 — azf~^){azf~^)^ 

and impose the symmetry constraint that Xij = Xj^i for i > j. Then the probability that a configuration 
[xij] gives a tableau of shape fi under the RSK correspondence is equal to 


c„(z, a, t)2:^i=i (t^^—t^'), 

l<2<j;<n 


(2.19) 


Cn{z,a,f) := z 


•= 7,-Ilj=i("-i)Q,-["/2] 


y-Y,f=iU-P{n-j) 


x{{i-azf-^) n (i-^v+^-^), 
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where hj := yij n — j and thus hi > ^2 > • • • > > 0. Furthermore in the scaled (Jacobi) limit 

t = e~^^^, z = a = hj/L = Xj, L—>■ oo, (2-21) 

when each lattice site (i,j) specifies a non-negative continuous exponential random variable with site 
dependent variance 

G [y,y + dy]) = (i + j - 2 + 2a)e“^^*+'’“^+^“^dy, i < j 

P'rixi^i G [y,y + dy]) = {i - 1 + a + i 

the probability (2.1!^ multiplied by L" tends to the PDF 


Cji (u, a \) . — 

where xi > X 2 > • ■ ■ > Xn > 0. 


r(a + oi + n) 1 



(2.22) 

n (e-^^-e-"‘), 
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r(2a + 2 z-l) 

2 = 1 
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After the change of variables and replacement of parameters 


Xn+i-j, a >—>■ (a + l)/2, ai i—>—Aj2, n 2n 


(2.25) 


we see that (|2.23D coincides with (|l.4[). 


Let us now consider the specialization (2.17) in (2.4). We must first give the form of (2.18) in the 
case that n i—> n 2 , A i—> /c, £(k) = ni with ni < n 2 so that Kni+i = ■ ■ ■ = k „2 = 0. Then with n* = ni 
and Tj := Kj + ni — j, manipulation of (^.1^) shows 






1 




n2 —ni —1 ni 

X (t; t)i 


(t, t)ri+n2 — ni 
(t]t)r, 


n 


(2.26) 


where the first product in the second line must be replaced by unity if n 2 = ni, ui + 1. Substituting this 
result, and (2.18) with n ni, n* = ni, A i—> /r in (2.4) we deduce an interpretation of the PDF ( |1.2| ) in 
the context of a last passage percolation model. 

Proposition 2. Let n 2 > ni. On each site of the ni x (n 2 + 1) square lattice specify a non-negative 
integer Xij according to the probability distribution 


Pr(cc,j=fc) = (1 - z¥+^-2)(z¥+^-2)^ j^n2 + l 

Pr(xi,r!,2+1 = ^) = (I - azf~^)(azt''~^)^. 


(2.27) 


The joint probability that a configuration [xi^j] gives, under the RSK correspondence, a pair of tableaux 
of shape pL, one of content ni and the other of content n 2 + 1, and that the subconfiguration [xij] i=i,...,ni 

j = l,...,n 2 

gives a pair of tableaux of shape k, one of content ui and the other of content n 2 , is non-zero if and only 

if 


hi > n > h2 > r2 > ■ ■ ■ > hni > rni > 0, 


(2.28) 


where hj := pij ni — j and rj := kj + ni — j. Furthermore the joint probability then has the explicit 
form 




(2.29) 






722—ni —1 ^1 ^2 

X n n 11 ( 1 - 


1=1 


i=ij=i 


i=l 


(2.30) 


In the Jacobi limit (2.21) (with the additional scaled quantity rj/L =: yj), (2.2t>) multiplied by tends 
to the PDF 


fcn„„.(a,ai)n(l-e-"^r”"^ H ( 

i—1 1<2<J<72i 
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(2.31) 
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^^1,^2 (^J^l) ■ — 


(.n?ll^im?li^l'-) r(a + ai) r( 2 a + z-l) 


n 


where it is required that 


xi > yi> X 2 > y 2 > ■ ■ ■ > Xm > Vm > 0. 

Analogous to ( 2.25| ), after the change of variables and replacement of parameters 

e 1 -^ X 2 n+i- 2 j 1 e X 2 n+ 2 - 2 j 1 a !—>■ {a + l)/2, oi 1 -^ —Al2 


(2.32) 

(2.33) 

(2.34) 


we see that with n 2 = ni = n ( ^.3l| ) coincides with ([l.Sl) . 

To specialize {2A) modified by the replacement ( ^.6|) according to ( 2.17| ), we note that for this to be 
non-zero (H) gives we require £{iu) < n 2 + I- Making the replacements k ^ y, n 2 ^ ni, ni ^ n 2 + 1, 
rj I—> hj := Kj + n 2 + ^ — j in (2.26|) shows 






n\—n2 —2 


n2 + l 


n (t;,), n n (c.-i-). 


Z=1 


(2.35) 




Substituting this result, and (2.18) with n i—> U 2 , n* = n 2 , A k, hj i—> rj := Kj + n 2 — j we deduce the 
analogue of Proposition || in the case ni > 712 . 

Proposition 3. Let ni > n 2 - On each site of the ni x (112 + 1) square lattice specify a non-negative 
integer Xij according to the probability distribution ( 2.2'^ . The joint probability that a configuration [xij] 
gives, under the RSK correspondence, a pair of tableaux of shape y, one of content ni and the other of 
content n 2 -f 1, and that the subconfiguration gives a pair of tableaux of shape k, one of 

content ni and the other of content n 2 , is non-zero if and only if 


hi > ri > h2 > r2 > ■ ■ ■ > hn^ > rn2 > h^^+i ^ 


> 0 , 


(2.36) 


where hj := p,j -\- n 2 1 — j and rj := Kj-\-n 2 -\-\—j. Furthermore the joint probability has the explicit 
form 


n2-ei 




X Yi 


(t',t) hi-\-ni — (n2-\-l) 
(t',t)h, 




(2.37) 


Kni,n 2 {a,z,t) := z 2 Ej=i ("2 j 
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1=1 
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(2.38) 
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In the Jacobi limit, ( 2.37 ) multiplied by i 2 n 2 +i tg ^/jg pjjp 


n2 + l 


i=l l<'i<_7<n2 + l 

Xe““5Z"=i(a:j+yj)g-ai T."j7Jxj-yj) ^-(a+ai)xn2 + l 
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where it is required that 
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711—712 — 2 
1 = 1 ‘ 


r(c 


ai 


(nr=r'^!)(nr=r'^!) r(a+ai) 


rii) r(2a + i + 712 — 1) 


r(2a + *- 1) 


(2.40) 


Xi > yi > X 2 > y 2 > ■ ■ ■ > Xn 2 > yn 2 > Xn 2 + 1 > 0. 


(2.41) 

gives the natural 


With ni = n 2 + 1, making the change of variables and replacements ( 2.34 ) in 
generalization of 0 to the case of an odd number of coordinates. 

In the case ni = n 2 =■ n, there is yet another combinatorial interpretation of the joint probability 
), which relates to a particular model (model (v)) introduced in 


Thus consider the 2n x 2n 
square lattice of sites (i, j), I < i,j < 2n. For the triangular shaped region specified by 1 < i < 2n — 1, 
i > j, i < 2n + 1 — j associate with each lattice site a non-negative integer Xij chosen according to the 
probability distributions 


Pi{xij=k) = Pr{x,^2n+i-j = k) = ^)'', 

Pr(a;i,i = fc) = {1 — azf~^){azf~^)^, 

Pl'(^7,2n+l —1—) — 9- 


i,j <n {i y^j) 


With the Xij in this region thus chosen, specify Xij at the remaining lattice sites in the square by the 
symmetry requirements 

Xj,i — ^i,j^ X2n+l — i,2n+l—j — ■ (2.42) 

The first of the symmetries in ( 2.42| ) implies that under the RSK correspondence the integer matrix maps 
to a single semi-standard tableau fi (of content 2n), while the second symmetry implies that each row is of 
even length. At a more sophisticated level, the resulting tableau is constrained to be self-dual (invariant 
under Schiitzenberger involution). Although we don’t present the details, using ideas from Q, from this 
one can show that with hj = p, 2 j-ij 2 + n — j and = y. 2 j 12 + n — j the probability that [xij] maps to 
the semi-standard tableau y is given by ( ^.29 ) with ni = n 2 = n. 

We remark at this point that the special case a = 1, A = 0 of the Jacobi parameter dependent PDF 
( |1.5D occurs as various probabilities in the work of Ciucu , and Krattenthaler on perfect matchings 
(tilings) on the Aztec lattice with removed sites. 

2.2 Laguerre limit 

It has already been remarked that after writing Xj i-^- \{xj 4-1), (j = l,...,2n), then scaling the 
variables and parameters according to ( |1.6| ) and taking the limit L —> oo, the parameter dependent 
Jacobi ensembles 0 and (0) reduce to the parameter dependent Laguerre ensembles ( |l.l| ) and (0 
respectively. As first noticed by Johansson |^, Laguerre ensembles can be obtained directly from the 


Schur measure (2.3) by first setting all the variables equal (and thus choosing t = 1 in (2.17)), then 
scaling the remaining parameters and variables as in ( ^.21 ). We thus obtain the following interpretation 
of the parameter dependent Laguerre ensembles. 
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Proposition 4. First, on each site {i,j) of the n x n square lattice, specify a continuous exponential 
random variable 


&[y,y + dy]) = 2ae dy, i<j 
Pr{xi,i G [y,y + dy]) = {a + dy 


and impose the symmetry constraint that Xij = Xj^i for i > j. Then the probability density that a 
configuration [xi^j] gives, under the continuous RSK correspondence of Appendix A, a non-intersecting 
path configuration with maximum displacement Xi at level-l, is given by (see also i/; 








{Xi Xj ) 


(2.43) 


where xi > X 2 > • ■ ■ > Xn > 0. For 2n this is equivalent to (1.1). 

Second, on each site (i,j) of the ni x (712 + 1) square lattice specify a continuous exponential random 
variable 


Pr(xij G [y,y + dy]) = 2ae dy, j ^ 712 + 1 
Pr(a;i„ 2 +i G [y,y + dy]) = (a + dy. 


(2.44) 


Then for 712 > ni the joint probability density that a configuration [xij] gives, under the continuous 
RSK correspondence of Appendix A, a non-intersecting path configuration with maximum displacement 
xi at level-l, and that the subconfiguration i=i,...,ni gives a non-intersecting path configuration with 


maximum displacement yi at level-l is non-zero if and only if the interlacing condition (2.3S) holds, when 
it has the explicit form 


nrir'/!nr=r'^! 

X n - Vj) 

l<2<j'<ni 


n 


no —ni 

2/i 


(2.45) 


In the case ni = 712 this is equivalent to the PDF (l.<). For < ni, the same joint probability density 


Vi, 


n 


2/z 


(2.46) 


is non-zero if and only if the interlacing condition ( 2.4-1 ) holds, when it has the explicit form 

X n n ~ yZ 

2.3 Gaussian limit 

It was pointed out by Baryshnikov upon interpreting a result of Glynn and Whitt |^, that for Xi^j 
i.i.d. random variables with finite variance, the quantity L{ni, n 2 ) specified by {2^) has a universal scaled 
form in the limit Tii ^ 00 , independent of the details of the distribution. This universal form is the PDF 
for the distribution of the largest eigenvalue in the GUE of 772 x 772 random complex Hermitian matrices, 
which have the joint eigenvalue probability density 


1 ^ 

^n. 


/=1 


(xk Xj) . 

I<j<k<n2 


(2.47) 
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It follows that with a = oi in ( ^.44 ) so as to obtain i.i.d. random variables, we can expect to obtain 
a Gaussian type ensemble by taking the scaled ni —> oo limit in the joint probability (2.46|). To see 


that this occurs requires nothing more than the classical transition between the Laguerre and Gaussian 
weights, 


a;i—»l+a;'y^2/c 


allowing us to derive the following result. 


Proposition 5. In {2.4-t) taking the Gaussian limit by setting 

a = ai, 2a = ni - (n 2 + 1) =: c, Xi^l+ Xi^/^Jc, yi\f^c, ni oo. 


gives the PDF 

2 n 2(«2 + l)/2 , 

^(n 2 +i )/2 n ® n ( 2 ;, - X,) n 

i—1 l<2<j'<n2 + l 

where it is required that 

00 > Xi > yi > X 2 > y 2 > ■ ■ ■ > > Xn^+i > -OO. 


(2.48) 


2.4 Limit to a biorthogonal Jacobi ensemble 


In the joint probability ( p.4D let us generalize the specialization (2.17) so that it involves two distinct sets 
of t variables and two distinct z variables, and thus choose 


(ai,...,a„J = (zi,ziti,ziq,...,zif”^ ), (61,..., = (z2,2:2^2, ^2*2: ^2^2^ )• 


(2.49) 


Using the Schur function evaluation formulas (|2.18| ) and ( 2.2(]| ) we can readily write down the generaliza¬ 
tion of (2.29|) and (2.37). Furthermore, the Jacobi limit of these generalizations can be computed. Let 


us make note of the explicit form in the case of ( 2.29 ). 

Proposition 6. Consider the generalization of l 2.2t\ ) obtained by specializing \2.f^ ) by ( 2.4ti J- Let 
zi = e~°-l^,Z2 = e~°‘/^Gi = /i^/L = Xj^rjfL = y^. 

Then as L 00 , this probability multiplied by Ip'^^ tends to the PDF 


2=1 




(2.50) 


^ni ,n2 (®’> *") ■ 


nr=i 1 ^/^! p(-Q, _|_ ^ _ X)+712) 


r(a + ai) 


n 

i=i 


r(a -f a -I- j — 1 ) 


(2.51) 


The Jacobi limit of the Schur function identity ( 2.15| ) tells us that if we integrate (2.50) over xi,... ,Xr, 
we obtain the ai —> 00 scaled limit (scaled by aj"”'^) of the same PDF, and thus 




-(o-l-a) Ujii Vj 


na 


— cyi\n2 —n\ 


2 = 1 


n 


-<^Vj _ _ p-Vi 


e-y^). 


(2.52) 


l< 2 <j'<ni 
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In the case rii = n 2 = n this same PDF was derived in ( A.10 ) of Appendix A from the continuous version 
of the RSK correspondence, giving the PDF for the event that the ni x n 2 lattice of non-negative random 
variables [xij] distributed according to 

Pr{x^j = y) = + -l) + a + 

gives a polynuclear growth model with height of the level-Z path yi. After the change of variables and 
replacement of parameters e~'^^ yn+i-j, (a -I- a) (a -I- 1) we obtain the PDF 

ni 

n - Vi) (2-53) 


2=1 


l<2<j?<ni 


where 1 > ?/i >•••>?/„> 0. For general a > —1, c > 0 and with m = n 2 the fc-point distribution 
corresponding to (2.53) has been computed by Borodin |M. As the method required to accomplish this 


task made use of ideas from the theory of biorthogonal systems, it was referred to as the biorthogonal 


Jacobi ensemble. The PDF (2.5C) represents a more general biorthogonal Jacobi ensemble, but the 
correlations for both sets of variables are yet to be computed. 


We remark that from ( 2.53 ) we can take a Laguerre and a Gaussian limit. The correlations for both 
cases were also computed in 0- 


2.5 Distribution functions 

Let us denote by i?(j; /; PDF) the probability that the interval I of the specified PDF (for this we will use 
the corresponding equation number) contains exactly j eigenvalues (for the sake of definiteness in termi¬ 
nology we will regard the PDFs as measures for eigenvalues). Similarly, let us denote by J;PDF) 

the same quantity except that only (•) =(e)ven labelled or (•) =(o)dd labelled eigenvalues are being 
observed. Let us suppose now that / = (s, oo) where s is inside the support of the PDF. Then as dis¬ 
cussed in [|^, knowledge of {i?(j; /; PDF)} is equivalent to knowledge of PDF)}. Furthermore 

p{k — l;s;ME) — the distribution function of the A:th eigenvalue from the right — is determined by 
{£;(};/; PDF)}. The quantities £;(};/; PDF), £;(')(};/; PDF), p{k - l;s;ME) for the PDFs (|^, (|^, 
(0, (0) are discussed in j^, as are the scaled limits of these quantities. Here we want to use knowl¬ 
edge of the so called hard edge scalings from 0 to identify scales associated to the large eigenvalues of 
the particular continuous RSK measures ( 2.23| ) and ( ^.3l| ). 


In (|1.4D, (O) we first change variables Xj i—> (I — X 2 n+i-j)/‘^ for consistency with 0. We know 


from 10 that then, with 2n = in (0 and n = N in ( |l.5[ ), the large eigenvalues have a well defined 
scaled limit obtained by setting 


1 

Xj = 1 - 

^ 2Af2’ 


A = 4N^a 


(we use a rather than a as used in |0 to avoid confusion with a as used in (2 .12| )) and taking N 
In particular 


lim 

JV —>-oo 

lim 

N —»-oo 

lim 

JV —»-oo 

lim 

N —>-oo 


A^4N^a 


E(-^(p;(l-: 




A=4N^S 


E{p- (0,s);SE‘'"’''^’“+1) 
E{p- (0,s);UE‘'“'^’“) 
E<-°\p; (0,s);OE^’“) 
E{p- (0,s); (OEU OE)“ 


OO. 

(2.54) 

(2.55) 

(2.56) 

(2.57) 
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Notice that ( 2.54 ), ( 2.55 ) are independent of the parameter a. We have already commented that the 
PDFs (2.23) and (2.31) (the latter with ni = n 2 = n) are related to 0 and ( |1.5D by a change of 
variables. Hence it follows that 


lim (p; (2 log — log oo); (|2.23|) 

lim (p; (2 log — log oo); (|2.3lD 
1^00 \ 4 


t=N,a^(a+l)/2 
ai = -2N^a 


l=n2 = N,a^(a + l)/2 
ai = -2N‘^a 


) = i?«(p;(0,s);(OE)^'“) 

^ = F;(:'(p;(0,s);(OEUOE)“’“)(2.58) 


Also of interest is the scaled form of the gap probability for the continuous RSK measure ( 2.52| ). Now 
results in imply 


lim E ( 

ni —»-oo 


/ q _ V f)P 

(?>; ( 0 . jVrrTTJ; (Si ) = 


4 ^1+1/c > 

\n2—niJ p\ d^P 
where is the integral operator supported on (0, s) with the kernel 




:= f; 


m—0 


m!r(a + bm) 


It follows from the relationship between ( ^.53 ) and (2.52) that 

lim E(p;((l + l/c)logAf-log^,oo);(|^) ) = -f>„ det(l - 

ni^oo V 4 -' "2-"i y p! ijfp 


«=i 


Notice that in the case c = 1 the scaled interval is the same as that in (2.5S). 


(2.59) 


(2.60) 


3 Interpolating ensembles from Macdonald polynomial theory 


The recurrence ( |2.14| ) satisfied by the Schur polynomials is a special case of a more general recurrence 
satisfied by the Macdonald polynomials, as is the marginal probability (|2.3|) and the evaluation formula 


( P.18|) . This then allows a generalization of the joint probability (^) to the Macdonald setting. We 
will see that the probability ( 2.19| ), which in the last passage percolation problem results from imposing 
the symmetry constraint Xi^j = Xj^i on the waiting times, is also a special case of the generalized joint 
probability. 


The recurrence (2.14) can be used to define the Schur polynomials. Likewise we can define the (monic) 
Macdonald polynomials Pk(6i, ..., bn^'^q, t) by the recurrence |^, pg. 348] 


(3.1) 


where with f{x) := {tx;q)ao/{qx;q)c 


ijf./K{q,t) := n 




(3.2) 

I I —j i,j " 1 ; ( d--'- r^J-r^ i.j " i 

The Schur polynomials are the special case q = t of the Macdonald polynomials (note that with q = t, 
f{x) = 1 and so ipfj,/K{q,q) = !)■ 


In the Macdonald theory, the generalization of the marginal probability (2.3) is 


Til Tl2 / \ 

n n (T ■ • ■, ; g, t)P^{bi, ..., g, t) 

*=ir=i 
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where 

and (•, 
is 


• j 0,ni j Qjt) — PfJ.) Pfj.iO'l ) • ■ ■ ) Orii j 0 


is a particular power sum inner product p8| , pg. 309]. Thus the natural generalization of (2.4) 

ni 712 + 1 , . 'i 

n n (T i ” Qm («1. ■ ■ ■: Q»1; 9 :(^1. ■ ■ ■. bn2 ; q, (g, ( 3 . 3 ) 

i=l j = l 

where the parts of /r and k are restricted by (2.5) in the case n 2 > ni, and (2.7) in the case ni > n 2 - 
Because of the identity dP), and the further Macdonald polynomial identity |^, special case of 6(a) 
pg. 352] 


n (a.b„.+i;g)oo ^ 

ira*o„2+i,gjoo 


(3.4) 


which generalizes ( 2.15 ), this furthermore has the interpretation as a joint probability density function 
on tableaux of shape ^ and shape k when /t/k is a horizontal strip. 

The Macdonald polynomials exhibit a generalization of the evaluation formula ( 2.18| ). Let := 
Pj denote the power sum of degree r. Define a homomorphism 


^u,t{Pr) — 


l-u’^ 

1 - r 


Then for any symmetric function / analytic in a:i,..., a;„ it is easy to see that 

£t".t(/) =/(l,t, ...,r"^), 


(3.5) 


(3.6) 


so the evaluation of Su,t{Pf_i) includes the generalization of (p.l8|) for the Macdonald polynomials. From 
pg. 343] we have with t = and £(A) < n 


!<“<„ “ (t;g)A.+fe(n-z) 

eu,t{Qx) = n (g^--^^-+V--^g)fcn /^( 

f^\{q-,q)x^+ki77-7) 


(3.7) 

(3.8) 


Note that for n fixed these quantities are strictly positive for u small enough, but may become zero or 
negative if n is unrestricted. In fact a special choice of u in each case shows that the action of £u,t may 
annihilate P\ or Qx for all A of length greater than a fixed size. Thus we see that 


et",t(^’A) = 0, egtn-i_t(QA) = 0, for £(A) > n. 


(3.9) 


while these same quantities are strictly positive for £(A) < n (in the first case corresponding to the 
evaluation (^)). 

In (|3.3|) we write &„ 2 +i Q;, replace by Ziai {i = I,..., ni), bj by Z 2 bj {j = 1,, 712 ) and take 

the limit ni,n 2 —>■ 00 . We would now like to apply the homomorphism t to the functions of { 0 ^} and 
the homomorphism to the functions of {bj}. For the Macdonald polynomials, after factoring out the 
z-dependence using homogeneity, this is done using the formulas (^) and (^). For the infinite products 
we recall the formula ]^ pg. 310] 


{zlZ2a^bj;q) 

00 


n {tZiZ2aibj;q)c 

1,7 = 1 


00 

n exp ( 


n—1 


(Z1Z2)" 1 - P 

n 1 — 9*^ 


■Pn(a)Pn(&)), 
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from which the action (^) immediately implies 




OO 

(n 


{ziZ2aibj\q) O 
{tziZ2aibj-,q) 

C 


OO 

-) = ne^p(- 


(Z 1 Z 2 )" (l-u")(l-w") - 
n (1 — — t") / 


n 


{ziZ2t^] q)oo{uwziZ2t'P\q)c 
{uziZ2tP] q)oo{wziZ2tP; q)c 


Aai} fr {zia i^n2 + l 5 ^)oo 
£u,t H (i2;iai6„2 +i;(3 ')oo 

Consequently the image of the joint probability 


p=0 

(zi6„^+i;g) 

OC 

{uzibn2+i]q) C 

is given by 


(MZia;(7)oo {uziZ2tP; q)ao{wziZ2tP; qj^o 


(3.10) 


For general u and w (3.1C) is not itself a meaningful joint probability on partitions /r, k because when 
the number of parts becomes large enough it will become negative. However, according to (|0) and ( ^.8[ ) 
for the special choice u = or w = this does not happen but rather ( 3.10| ) vanishes when £(^) > n. 
Thus we are naturally led to two distinct joint probabilities on partitions /i, k with /t/k a horizontal strip, 


Pre(/r, k) := 


[zia]q)c 


n 


{ziZ2tP ; q)oo {wziZ2tP +'^; q)c 


for which 


and 


(t”zia;q)oo {ziZ2tP+'^\q)oo{wziZ2tP]q)o 
X £t-,t (Q/x )e«,,i (P/, )'!/'m/«; (l- 0 ^ 1 ^'4'^'" I'‘I 

fbl > Kl > fJ,2 > K2 > ■ ■ ■ > fJ-n > Kn > 0, 


(3.11) 

(3.12) 


Pro(^, At) := 


{zia-,q)o 


n 


{ziZ2tP;q)oo{uziZ2tP+'^ ^■,q)c 


{uZia;q)oo ^p^{uziZ2tP]q)oo{ziZ2tP+”^ ^<?)oc 


(3.13) 


for which 


A^l > Atl > /r2 > A«2 > • • • > Mn-l > At„_l > > 0. 


(3.14) 

A short calculation shows that setting w = zi = Z 2 = z, t = q in Pre(/r, k) gives ( 2.29D , 

while setting u = , zi = Z 2 = z,t = g in Pro(Ai,At) gives (|^). Furthermore (|^, (|^ 

in the case t = q^ reclaim ( 2.19|) . Thus straight forward simplification gives the following result. 

Proposition 7. Let w = q~^t'^, zi = z, Z 2 = qz, and t = q^, and write h 2 j-i '■= tij + 2n — (2j — 1), 
h 2 j ■= Kj + 2n — 2j. Then (3.11) reduces to (2.It) with n 2n, t q in the latter. Similarly, let 
u = qt'^~^, Zi = z, Z 2 = qz and t = q^, and write h 2 j-i := ^.j + (2n— 1) — (2j — 1), h 2 j ■= Kj + (2n— l) — 2j. 
Then ( 3.IS ) reduces to \2.1!^ with 2n — 1, t ^ q in the latter. 

The probabilities Pre(/t, k) and Pro(/r, k) exhibit a special property with respect to summation over 
/r. Thus it follows from (3.4) that 


y]pre(/i,K) = jq 


{ziZ2tP; q)oo{wziZ2tP~^'^; q)c 
pLf iziZ2tP+'^-,q)ooiwziZ2tP;q)c 


,t (,Q K)^W,t )(yZlZ2 


iIkI 


V-Ti A ^ iziZ2tP-,q)^{uziZ2tP+^ ^■,q)oo ,|«| 

Z^Pr°(/^,«) = 11 euAQ.)et^-^,tiP.)iziZ2y ' (3.15) 


P—0 


{uziZ2tP-,q)ao{ziZ2tP+'^ ^■,q)c 
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A special case of Pro(/i, k) also exhibits a special property with respect to summation over k. To see 
this we first note from (|3.6|) that 


and so for the k dependent terms in ( |3.13| ) with a = 1, Z 2 = t we have 

K K. 

= P^(l, r-i) = (3.16) 

where the second equality follows from (|3.lD and the fact that P^ is a symmetric function. Thus 


^Pro(/r, k) 


{zi;q)oo {zitP+^-,q)oo{uzitP+^;q)oo _ 

(uzi;g)ooj-ijw^itP+i;g)oo(^itP+";g)oo ' 


(3.17) 


To use (3.16) in the case of Pre(/i, k) we must set w = This in turn implies k„ = 0, so we see 


that no new identity results, but rather we reclaim the special case u = t" of (3.17). As made explicit 


in Appendix B, the identity (3.17) can be recognized as being equivalent to a special case of a g-integral 
due to Evans |l^, and also as the v = % case of Okounkov’s q-integral representation of the Macdonald 
polynomial ill- The structure afforded by (|3.17| ) suggests a simplified derivation of the latter which 
is given in Appendix B. We remark too that ( ^.17 ) can be considered as a particular q,t generalization 
of a class of measures on partitions known as z-measures § . 

Because the probabilities ( 2.29| ), ( 2.37 ) and ( 2.19 ) can all be derived from ( 3.1C ), they all exhibit 
the special property ( 3.15| ). Of course in the cases of ( 2.29| ), (2.37) this identity is immediate from their 
interpretation as joint probabilities for tableaux of shape fj, and tableaux of shape k with /x/k a horizontal 
strip. But the probability (p.19) has no such interpretation, and the identity implied by (3.15), 


even 


l< 2 <_j<n 

Y[ 




Z^j = l 


a=0 


^2j —1—(j —1,..., [n/2]) 1 
hji=Q (n odd) 


(3.18) 


(here the notation even( ) denotes the distribution of the even labelled coordinates h 2 ,h 4 ,..., /i 2 [n/ 2 ])) 
telling us that there is no dependence on a after summing out the odd labelled coordinates, cannot easily 
be anticipated. 


3.1 Pto{p,k) and Pre(p, k) in the Jacobi limit 

Consider (|3.13|) with zi = z, Z 2 = tq~^z, u = t = q^. The Jacobi limit is obtained by setting 


— p-afL 


z = e "'“,z = e 

2n-l 


— ajL 


a = e “i /^,9 = e , fij / L = Xj, Kj / L = yj, 


multiplying (3.13) by ^ and taking the limit L ^ oo. This gives the PDF 

n 

Cn{a, a, ai,f3, A:)e““^^=i 


2-“Er=i yig-“i(E?=i 2:i-E”=i Vi) P(1 — e“^*)^ 


X n 

l<2<j<n 


n I® 


2=1 

n n—1 

-Vi _ p-vi \ TT TT \p-^i - p-vi\k-i 


! _ ^ r(a + tti + /3 + 1 + A:(n — 1 )) -p-|- 

C„(o,a, ai,/3, fc) :=- ^ ^-n/a , - H 


iniT' 

2=1 i=l 
n—1 


(3.19) 


r(a + ai)r(/3 + 1 ) 




r(a + a + /3 + k(ri — 1 + f)) 
T{ki)r{P + ki + l)r(a + a — 1 + ik) 
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Replacing a,abya + l,a+l and changing variables e 

n 

Cn{a + l,a + 1,01,/3, A:) (1 - Xi)^ y[ 


a;„+i_i, e 1 -^ yn-i, this reads 


n—1 


a—ai 

3 


n 


i=i 

n n— 1 


n \y3-y^\W_W_\^3-Vi 


k-1 


l< 2 <j<n 




i=l 


and we require the analogue of the interlacing condition (3.14), 

I > Xi > yi > X 2 > y 2 > ■ ■ ■ > y-n-l > Xn> 0- 


(3.20) 


(3.21) 


A natural generalization of ( 3.2C ) is to include a factor 0^=1 i^~yj)^- To compute the corresponding 
normalization, we note that with R denoting the region ( |3.21 ), the Jacobi limit of the second identity in 
( tl.lS ) tells us that 


dxi 


IR 


■dxnWxf(l-XiY \Xt-Xj\Y\Y\\xj - y^ 

i—1 i—1 j—l 


k-1 


r(i + Q)r(i + /3)(r(fc))"-i 

r(2 cx -\- j3 -\- {ji — 1)^) 


n—1 

i—1 


ly^-yjl 


2k-l 


(3.22) 


This is an integration formula due to Anderson |^] . Furthermore, we have the well known Selberg integral 
evaluation 


Jo Jo 


|2A 


i<j<fc<Ar 


N-1 

n 

3=0 


r(Ai +1 + j A)r(A 2 +1 + jA)r(i + {j + i)A) 
r(Ai + X 2 + 2 + {N + j- l)A)r(l + A) 


=:5iv(Ai,A2,A). (3.23) 


It follows from ( ^.22 ) and (3.23) that 

jCn,™-!)/ N r(2 + a + P + {n- l)k) 

O ’ y / ’ -pi / -| I \ pi / -I I /Q\ / T~' / 7 \ \ n 


r(l + a)r(l + /3)(r(A:))—1 Sn-i{a + ai + k, P + Pi + 


n—1 


n n—1 


X JI \xj-Xi\W_y'^^{l-y,f^ JJ \y3-y^\W_W_\x3-y^\^ ^ 

l<i<j<n i—1 1 i—1 j—1 


(3.24) 


is a correctly normalized joint PDF. With a = a + ai, ai = a — ai, /3i = 0 it coincides with (3.20). 


Key properties of ( 3.24 ) are 

f dxi--- dxn {x, y) = - - 

JR <t>n—1 


1 


n—1 

X n n \y^ - yj\"" =■ 

i—1 1 


(3.25) 


which follows immediately from ( 3.22| ), and 

[ dyi ■ --dyn-i Ji”’”"^^(x, y) 
JR 


ai=/3i=0 


1 " 

-ir^Yix^l - Xif JJ \xi- Xjl^’^ =: JIP^'-\x). 


Sn{a,P,k) J 


(3.26) 


Z=1 


l<i<j'<n 
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which can be deduced from the Jacobi limit of ( 3.17 ). Like (|3.22| ) (and thus ( 3.25| )), ( 3.26 ) is an integration 
formula due to Anderson [01. 


Consider now ( 3.11 ), and put Zi = z, Z 2 = tq~^z, w = t = q^. The Jacobi limit is obtained by 

scaling the parameters according to ( 3.19 ), multiplying ( 3.11 ) by and taking the limit L ^ oo. One 
thus obtains the PDF 

n 

A:„(a, a, ai,P, fc)e-“^?=i 


X n I 

l<z<j'<n 


Kn{a,a,ai,P,k) := 

By changing variables i—> 

the PDF 


i=l 

^^-^3 — e-V3\k-^ 

*,1 = 1 

F(a + Cl + knj) 1 


n - e 


n+h 


n 


F(a + a - 1 + /c(j + 1 + n)) 


r(a + ai) +i)) jJi T{a + a - I + kj) 

Xn+i-i, e~y' I—> Un-i and replacing a, a by a + 1, a + 1 we obtain from this 


A:„(a+l,a+l,ai,/3,A:) (1 - \xj - Xi\\yj - y^\ \xj - y^\^ ^ (3.27) 


l<2<j<n 




and we require the analogue of the interlacing condition (3.14), 

I > yi > Xi > y2 > X2 > ■ ■ ■ > Vn > Xn > 


(3.28) 


Unlike the situation with ( 3.20 ), it is not a natural generalization of ( 3.27 ) to include an extra factor 
(here nr=i(l — XiY). Only with this extra factor absent can we compute the integral of the x's according 
to the Jacobi limit of the first identity in ( 3.15 ). Before stating this result, let us first rename some of 
the parameters and manipulate the normalization so that ( 3.27 ) reads 

jO,«)x N r(l + a + nfc) _1_ T\ x^v^^Hl - 

" T{l + a){T{k)r Sn{a + ai+k,Pi,k)l^ ^ 


X la;^- - x,\\yj - y*| H 


k-1 


(3.29) 


In terms of this quantity the Jacobi limit of the first identity in (3.15) reads 

1 


[ dxi - ■ ■ dxn Ji^''^\x, y) = 
JR 


Sn{a + ai + k,(3i,k) 


i—1 

where R denotes the region (^.28). It follows from this that we also have 


(3.30) 


dyi---dynJ^"’"'\x,y) 


1 


ai=o Snia,Pi+k,k) 


Y[xni-x.f^+^ n \x^-x,\^'^=: 4 -’-\x), 

i=l 


(3.31) 
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4 Random matrix interpretation of the Anderson density 


With Jo"’” ^\x,y), J^’'\x) and ji"’” ^\y) defined by ( f3.24 ), ( 3.26| ) and ( 3.25D respectively, we can 
construct the conditional PDF’s 


r(nfc) ni<j<j<„-i(2/i “ %) 
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where the a;’s and y’s are interlaced according to (3.21). The conditional PDF (}4.lD can be recognized as 


the special case si = S 2 = • • • = = A: of the conditional density function 
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appearing in the work of Anderson Q on the Selberg integral. (Note that the ^-generalization of this 
same density appears in Evans’ q-integral ([B.5|); also we refer to for a further integration method to 


verify that the normalization is correct.) Here we will show that (O), with the Si non-negative integers 
or half integers, and thus according to © the conditional PDF associated to the case a = oi = 0 of the 
interpolating ensemble ( 3.20| ), can be derived from a random matrix problem. 

The random matrix problem relates to the corank 1 random projection of a fixed matrix. Specifically, 
we seek the eigenvalue PDF of 

M := HAH, n := 1 - xx^ (4.4) 

where A is a real symmetric, or complex Hermitian, fixed matrix, T is a real, or complex, normalized 
Gaussian column vector of the same number of rows as A and 1 denotes the identity matrix. The 
eigenvalue PDF depends only on the eigenvalues of A, which we take to be oi > 02 > • • • > a„ with 
multiplicities mi, m 2 , ■ ■ ■, mn- 


All but n eigenvalues of (4.4) must coincide with the eigenvalues Oi of A and must occur in M with 
multiplicity — 1. For the latter result we make use of the following formula for the characteristic 
polynomial of M. 


Lemma 1. We have 


det(M — Al) = —Adet(A — Al)Tr^(A — Al) ^xx^^. 


(4.5) 


Proof. Simple manipulation using {^A) shows 

det(M — Al) = det(A — Al) det(l -I- (A — Al)~^(—Aif:?’*’ — xx^A -I- xx^ Axx^)). 

The matrix in the second determinant is of the form 1-l-P where Y has rank 1, and in such a circumstance 
we have det(l +Y) = 1 -I- TrP. Using this fact, then further manipulation using Tr(xx^) = 1 (which in 
turn follows from the assumption that x is normalized) gives (4.5). □ 


Now (4A) shows there is an eigenvalue A = 0, and that the remaining eigenvalues satisfy 


n n (j) 

- A)’"' : = 0, 


Z = 1 


ai - X 




0) 


(4.6) 


19 



























where the denote the diagonal elements of xx^. This shows immediately that M has eigenvalues 
with multiplicities 771^ — 1, and the remaining n — 1 eigenvalues given by the zeros of the random rational 
function 


R{\) := 


Wi 


i=l 


a* - A 


(4.7) 


The fact that the Wi are positive (being equal to sums of squares) implies that the roots of R{X) are all 
real (as must be since M is Hermitian) and further have the interlacing property 


Oi > Ai > 02 > A 2 > • • • > An— 1 > Or! 


(4.8) 


(c.f. (^.2lD). 

We would like to compute the distribution of the roots of A for given oi,..., o„ and the {wi} random. 
This depends crucially on the precise distribution of the {wi}. Now Wi has the form 


Wi — Xij{Xi + • • • + Xn) 


(4.9) 


where Xi consists of /3mi (/3 = 1 for x real and /3 = 2 for a? complex) independent real Gaussians with 
mean zero and standard deviation cr, and thus has the gamma distribution r(si,2(T), Si := l3mi/2. It 
follows that the PDF for {wi,W 2 , •. •, Wn-i', Wn) is equal to the Dirichlet distribution 


r(si + • • • + Sn) 7 - 1 - Si-1 

r(.i)---r(s„) n-. 


n—1 

Wn ■= 1— ^ Wj, Wj > 0. 


(4.10) 


The working in Anderson’s paper Q| shows us that the distribution of the roots Xi of R{X) when the 
{wi} are distributed according to (4.10) is given by (4.3) with yi = Xi, Xi = ai. As a consequence we can 
specify the sought eigenvalue distribution. 


Corollary 1. The eigenvalues Xi,... ,Xn-i of M in & differing from the eigenvalues of A and from 
0 have the PDF with yi = Xi (i = 1,... ,n — 1), Xi = Oi (i = 1,... ,n) and Si = j3mil2 (/3 = 1 for 
X real and (3 = 2 for x complex). 


In the case that all eigenvalues of A are distinct (or doubly degenerate in the case of x complex), 
the eigenvalues of M can be interpreted as so called radial Gelfand-Tzetlin coordinates introduced by 
Guhr and Kohler [^. Further, as shown in Appendix C, this observation and Corollary can be used 
to rederive a recursion formula obtained in for certain matrix Bessel functions. 


4.1 Construction of interpolating Jacobi ensembles 

We can make use of Corollary to determine explicit random matrices with eigenvalue PDFs which 
realize (Fd) and ([l.Sf ). Consider first (H- Essential to our construction are random matrices with 


a doubly degenerate spectrum which have an eigenvalue PDF of the form (1.9). The required random 
matrices are known from Q] (see also [^). Thus consider a member S of the circular ensemble CSE(„._|_„) 
(for the definition and construction of such matrices — in which each element is itself the 2x2 matrix 
representation of a real quaternion — see e.g. @). Decompose S as 


S = 


'^n* Xn 

^nxn* 


^n* xn 


(4.11) 


where n* > n and the subscript on the blocks tells us their dimension (with the already mentioned 
qualihcation that each element is a 2 x 2 matrix). Then we have from 0,0 that the random matrix tt^ 
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has eigenvalue PDF (O) with a = 2(n* — n). Moreover, if we append to t an extra no rows of zeros and 
denote this i say, then iP has a zero eigenvalue of multiplicity ng, and n eigenvalues with PDF ( |1.9| ) of 
multiplicity 2. Substituting iP for A in (|4.4| ) we can make use of Corollary to deduce that the PDF 
for the non-zero eigenvalues of M realize 0 )- 

Theorem 1. With t specified above and x a normalized complex Gaussian column vector of2n + no rows, 
the non-zero eigenvalues of the random matrix 


M = m?n, 


T\ = l — XX 


have the PDF 0 with 


a = 2{n* — n), A = 2{n* — n) — 2no 1. 


(4.12) 


Proof. Let the non-zero eigenvalues of iP be denoted Ci, 02 ,..., a„. We know they are doubly degen¬ 
erate and have distribution 0) with a therein given by ( [4.12[ ), and furthermore we know that iP has a 
zero eigenvalue of multiplicity uq. According to ( [4.6| ) and the matrix M then has a zero eigenvalue 
also of multiplicity no, eigenvalues oi,..., with multiplicity 1 and eigenvalues Ai,..., A„ such that 


fll > Ai > 02 > ■ ■ ■ > An > 0. 

From Corollary the conditional PDF of Ai,..., A„ given oi,..., a„ is proportional to 
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(4.14) 


But the PDF of ai,...,an is given by ( |l.9| ) (with the x’s replaced by a’s and a given by ( [4.12D ), so 
forming the product with ([4.14) shows the eigenvalue PDF of M is proportional to 


n- 

2=1 


2(n*—n) — 1 / \'^0 ^ 


(~) n n 


After relabelling we recognize this as the PDF ( |1.4[ ). 

A significant feature of Theorem is that by construction 


odd(Af) = JSE,, 


6=0 


□ 


(4.15) 


where odd(M) refers to the distribution of the odd labelled eigenvalues of the random matrix M. Choosing 
any particular value of A allowed by (4.12) and multiplying both sides of ( l.lSj ) by nr=i 
that (jnsj) implies 

odd(JOE2„U=(c-i)/2) = JSE„| a=c+i (4.16) 


where c = 2no — 1. Identities of this type were classified in using functional properties of the 
PDFs. In fact it was found that (4.16) is one of only two identities relating every second eigenvalue in 
a matrix ensemble with orthogonal symmetry and an even number of eigenvalues, to a matrix ensemble 
with symplectic symmetry (for the other see ( |5.9| ) below). A challenge was issued to provide a matrix 
derivation of such results; by way of the above construction this challenge has been answered for the 
particular identity ( 4.16 ). 

Let us now seek a realization of 0 as an eigenvalue PDF. Guided by the above construction of 
random matrices with eigenvalue PDF (|l.4[), we first seek random matrices with a doubly degenerate 
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spectrum which have an eigenvalue PDF equal to ( 1.12| ) which is the A —oo limit of (|l.4| ). From 1,0 
we know that with S a (n* + n) x (n* + n) random unitary matrix decomposed as in ( |4.11| ), the random 
matrix has eigenvalue PDF (1.12|) with 


a = n — n, 


(4.17) 


although the eigenvalues are all distinct. To obtain a doubly degenerate spectrum with the same eigen¬ 
value PDF, we simply replace each complex element a; -I- iy of t by its 2 x 2 real matrix representation 


X y 

-y X 


(4.18) 


To this doubly degenerate spectrum with eigenvalue PDF (1.12) we want to add a zero eigenvalue of 


degeneracy uq. As noted below (4.11), this is achieved by simply appending no rows of zeros; let us 
denote the real representation of t so modified by i. The real symmetric matrix iP" then has a zero 


eigenvalue of multiplicity no and n eigenvalues with PDF (1.12) of multiplicity 2. We can now use 
Corollary to obtain the sought realization of ( |l.5| ). 

Theorem 2. With t specified above and x a normalized real Gaussian vector of2n+no rows, the non-zero 
eigenvalues of the random matrix 


M = mFn, 


n = 1 — XX 


have PDF O with 


a = n — n. 


A = n* — n — no-\-l. 


Proof. Following the reasoning of the proof of Theorem |^, the matrix M has a zero eigenvalue of 
multiplicity no, a distinct copy of the non-zero eigenvalues oi,..., a„ say of iP", and eigenvalues Ai,..., A„ 
satisfying the interlacing condition (4.13). Corollary ^ with n i—> n -I- 1, a„+i = 0, s„+i = no/2, Si = 1 
{i = 1,... ,n) gives that the conditional PDF of Ai,... A„ given oi,...,On, is proportional to 


-i/A*Wo/2-i 
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n 


The eigenvalue PDF of M now follows by multiplying this by the PDF of ai,..., o„ as given by (1.12) 


(with the x's replaced by a’s). After relabelling the coordinates the PDF ( |l.5[ ) results with the parameters 
as stated. □ 


Analogous to (4.15), by construction 


odd(M) = JUE, 


6=0 


(4.19) 


Since with A = 0 the eigenvalue PDF of M coincides with that of the matrix ensemble JOE„| a=(a-i )/2 U 

b=0 

JOE„| a=(a-i )/2 (recall sentence below (Oil)) we have a matrix theoretic understanding of the relation 


OddfjOE„|a = (a-l )/2 U JOE^I a = (a—l )/2 — JUE^ 

V 6=0 6=0 / 


6=0 


(4.20) 


4.2 A random three term recurrence for interpolating Jacobi ensembles 

In this section, inspired by the recent work |^, it will be shown that J^'^~^\x,y), specified by (3.24), 
and Je^’^\x,y), specified by ( 3.29|) , can be sampled from the zeros of a polynomial which in turn is 


22 

























specified using a random three term recurrence. Consider first (B- We begin by noting that ( [4.21) , 
like ( [4.1[ ) is intimately related to the Anderson density (4.3). Thus in the latter put n n + 1, relabel 
the x’s by y’s and the y’s by x’s, then set yi = 1, y-n+i = 0, relabel j/i+i hy yt {i = 1,... ,n — 1) and 
put Si = /3 + 1, s„+i = a + 1 and Si = k {i = 2,... ,n) to obtain ([4.2|) . As a consequence, Anderson’s 
result relating the random rational function (E3)’ with coefficients distributed according to the Dirichlet 
distribution (4. H), to ( [4.3D tells us we can similarly specify a random rational function related to ( [4.2| ). 

Corollary 2. Denote the Dirichlet distribution by Dn\si ,..., s^-i; s„]. Let {wq, ... ,Wn-i;Wn) be 

distributed aceording to + 1, (fc)"'“^;a + 1], where the notation (fc)"“^ denotes k repeated n — 1 

times. We have that the roots of the random rational function 


.^n+l(^) ■ — 
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Wn 
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h^-y^ 


(4.21) 


are distributed according to the PDF 0- 

Anderson’s result stated below ( j4.1(]| ) and Corollary ^ can be used to derive a random three term 
recurrence which specifies a polynomial, the zeros of which sample from the joint PDF ( 3.24 ), and also 
sample from the marginal PDF (3.26). We will specify the recurrence by first detailing how it leads to a 
polynomial with the sought properties in the low degree cases, before stating its general form. 

Step 1 Consider ( [4.2lD in the case n = 1 and write wq , wi . Let (wq^^ ; ) be distributed 

according to + 1; + 1]. Let denote the zero of (4.21) in this case and form the polynomial 


Ai(x) := X — 


( 1 ) 


It follows from Corollary ^ that A^^^ is distributed according to 


{x,y) 


{x) 




=: P(A«) 


(4.22) 


(4.23) 


which is itself the Dirichlet distribution D 2 [aA'> + 1; (3^) -\- 1], 
Step 2 Define Ai{x) by ( [4.22 ) and also define 

Aq{x) := 1. 


(4.24) 


Let be distributed according to + l,k;a^'^^ + 1] and construct the random 

quadratic polynomial 


A 2 (x) := w^\x — \)Ai{x) + Wq‘‘'’xAi(x) + wl'^^xlx — l)Ao(x). 


( 2 ), 


„( 2 ). 


(4.25) 


Dividing both sides by x(x — l)Ai(x), this reads 
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X - Af' 


x{x — l)Ai{x) X x—1 

^ o f2l f2') 

Because 1 ^^^e wf ' are positive, A 2 (x) must be monic with real roots and we write 

A2{x) = {x - X^i^){x - X^2^). 

It follows from Corollary^ that the conditional distribution of {aJ^\ given A^^^ has the form 
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(4.27) 


and so the joint density of is 

P(Afy(Al"\Af |Af^ 

If we set 

0,(1) ^ ^(2) +ai+k, /3(i) = /3(2) +/3i + k (4.28) 

we recognize this as the joint distribution function J^'^\x,y) with a = a^'^\j3 = According to 

C^“^Tv f2l ( 2 '\ 

(p.26), if we now set ai = /3i = 0 we can compute the marginal distribution of A) ^ A 2 ^ 


A®>AW>Af) 


dAWp(AW)P(A('\A(")|A«) 


ai=i3i=0 


= j(^-)(Af)*,A(^)*) 


« = a(2) ’ 

/3=/3(2) 


(4.29) 


where the use of the * on the right hand side indicates the parameters have been chosen so that (4.28) 
holds with ai = Pi = 0. This distribution is realized by the roots of ^ 2 ( 0 ;) when A$^^^ is not observed. 


To proceed further requires an extension of the above arguments. Let us consider 
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Ai{x) 

A*(x) 


= E- 

1=1 


Ui 


-a: 


( 0 *’ 


Ul + U2 = 1 


(4.30) 


and pose the question as to what distribution of (ui;u 2 ) is required so that the distribution of the zero 

fll (2')* f2^* 

on the right hand side has the same distribution as A) , with A) ' , A 2 distributed by the right hand 
side of ( 4.29 ), and is thus specified by ( 4.23 )? 

According to Anderson’s result stated below ( 4.10| ) we have that with (ui\U 2 ) distributed according 
to D 2 {k\k), the root of (4.30) has conditional distribution 


t(24)q(2)* ;v(2)*, >(1)n 
Jo , A2 , Ai ) 


ai=/3i=0 


Thus the corresponding marginal distribution of A^^^ is 


7 ( 2 ,1)^(2) \(2).\(1)'> 

Jo , A 2 , A]^ ) 


c. = ct(2) 

f3=l3(^') 


ai—f3i—0 




= J^1)(A«) 


c. = a(2)+fc 
y3=/3(2)+fc 


where use has been made of ( tl.25| ) and ( 4.28 ). This is indeed the same distribution as ( |4.23 ), provided 
we set ^( 1 ) = ^( 2 ) _|_ therein. Let us denote Apx) with the parameters so specialized 

by Af{x). 

As well as making use of ( 4.30 ) with (iti, U 2 ) distributed according to D 2 [k\ fc], we require some special 
properties of the Dirichlet and beta distributions. First we recall that the Dirichlet distribution D 2 \a] P] 
and the beta distribution B[a,P] are the same thing. We require the fact that if (wq, ■ ■ ■, Wn-i; Wn) is 
distributed according to Dn+i[ao ,..., Un-i', Qf„], then the marginal distribution of Wj (j = 0,..., n — I) 
is given by B[aj,J27=o i^j ^i] ^’^d the marginal distribution of Wj + Wk, (j k, j,k < n) is B[aj + 
i^j k require the property of the beta distribution (see e.g. [^, pg. 42]) 


B[a + b, c]B[a, b] = B[a, b + c] 


(4.31) 


where here — in an abuse of notation — the left hand side means the product of random variables from 
the respective distributions, and the right hand side tells us the distribution of the product. 
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Step 3 Analogous to the construction of ^ 2 ( 0 ;), we construct A 3 (a;) by the random three term recurrence 


A 3 (x) := w^\x — 1)A2 {x) + wl^'^xA^ix) + w^^'^x{x — l)Af (x) 


where {wl^\ ) is distributed according to I? 3 [/ 3 ^^^ + + 1 ], or equivalently 


Asix) 


„(3) 


x(x— 1)A2 (x) X 
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(3) 

X — 1 


.,(3) 


At{x) 

A*(x) ■ 


For Af (x)/A 2 (x) we substitute ( |4.30| ). Now the theory above ( [t.31| ) tells us that the marginal distribution 
of is i3[2fc, + 2], while the distribution of ui in ( 4.30 ) is B[k, k]. Applying ( 4.31 ) it follows 

that we can write 

, ( 3 ) Af (x) Ui 


w 


1 A*(x) 


= E' 




where ui has distribution B[k,a^^^ + + fc + 2 ] and mi + U 2 has distribution B[2k-,a^^^ + P^^^ + 2 ]. 

Consequently we have 


A 3 (a:) 


r,(3) 


r,(3) 


E- 


r,(3) 


X — 1 ^—' X — A 

X A2 


( 1 )* 


(4.32) 


x(x — 1 )A 2 (x) 

where w^\ ) has distribution Di[P^^'> + 1 , (kp ; + 1]. Arguing now as in the derivation 

of ( 4.26|) that 

A3(x) = (x - Ap^)(x - \i^^){x - 

where the conditional distribution of A^^^ A^^^} given {A^^^*, has the form 




(-,2) 


(y) 


— --^ 1^1 5^2 5^2 ’^2 )' 


The joint density of {A^^*, Ap\ A^^^ A^^^} is therefore 

j(2A(Af)*, A^'^*)P(Ap\ Af, A® |Ap^*, Af *) 


(4.33) 


which with 

q,(2) _ q,(3) _|_ _|_ ^^ pi2) _ p(3) ^ 

we recognize as the joint distribution function Jo^’‘^\x,y) with a = P = P^^\ As with (|4.29|), the 

(3) (3) (3) 

marginal distribution of A) , A 2 > Ag can be computed in the case ai = Pi = 0. Thus it follows from 
( ^.26 ) that 




dAf dAf f 2’-)(Af Af *)P(Af, Af, Af |Af *, Af *) 


ai=/3i=0 


= jf-)(Af*,Af*,Af*), 


(4.34) 


giving the distribution of the zeros of ^^(x) when Af *, Af * are not observed. 

Step 3 is representative of the general step n in generating the recurrence. Of course we now need 
inductive hypotheses relating to the roots of polynomials generated in earlier steps. In particular, we 
suppose that in step n — 2 a polynomial A„_ 2 (x) has been generated and the density of its roots is given 
by J^-’"“^Hy)lai=/ 3 i=o as specified by ( $.25[ ) with a = P = With the special choice of 

parameters + k,P^‘^~^^ = + fc we denote A„_ 2 (x) by A^_ 2 {x). At step n — 1 
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we require that a polynomial A'^_i{x) has been generated which has the density of its roots given by 
specified by ( 3.26| ) with a = , (3 = 

Step n We construct An{x) by the random three term recurrence 

An{x) =w^J"\x- l)^*_i(a;) + xAl_^{x) +w^^'’xix- l)Af_2{x) (4.35) 

where {w^\ w'f'\ ) is distributed according to D3+ 1, (n — + 1]. Arguing as in the 

derivation of (4.33) we see that with 

+ai+k, =/?(’")+/3i + fc (4.36) 

the joint distribution of the roots of An{x) and A"^_^ (y) is given by (x, y) with a = , /3 = , 

and that with ai = /3i = 0 the marginal distribution of the roots of An{x) is given by Jo”'’~\x) with 
a = a(’"), 0 = 

From a practical point of view, our objective is to sample from Jo"’" ^\x, y) and Jo^’~\x) for a fixed 
value of n and fixed parameters. To sample from Jo^'~\x) with a = oq, /? = /So we implement the above 
steps with 

Q,0) — pU) — j'jj^ _|_ (4.37) 

We see that in this situation Af{x) = A*{x) and so {A'f {x)}j=2,...,n is determined by the random 
recurrence 

A"^ (x) = W 2 \x — l)A^_-y{x) + Wq^ xA^_-^{x) + w^^x{x — l)A^_ 2 {x) (4.38) 

where {wq\w^\w 2 '^) is distributed according to D^[{n — j)k + /3o + 1, (j — l)fc; {n — j)k + q/q + !]■ 
The initial conditions for the recurrence are A^-^ix) = 0 and A^{x) = 1. The zeros of A^{x) then are 
distributed according to Jo"”^(a^)| “=“o. 

If our objective is to sample from ji"’" ^^(x, y) with a = ap, (3 = /3o, we again compute {A^ (a:)}j=o.....n-i 
this time replacing ag, /3o by oq + ai, /3o + /?i throughout. Let us write A"^{x) with these parameters as 
(cc). Because we now have = fc + op + ai, = k + Po + Pi we see that A^_^{x) = A^_^{x), 

so according to (4.35) the final step is to compute 

An{x) = w^\x - l)A*_-y{x) + w^^\a*_^{x) +w^^'’x{x- 1)A*_2{x) (4.39) 

where {w^\w^\w^'^) is distributed according to D^[j3Q + 1, (n — l)fc; oq + 1]. We then have that the 
zeros of (A„(a;), A^_^(y)) have the joint distribution Jo"’"~^^(a;, y). 

Let us now turn our attention to sampling from (a:, y) as specified by (3.29). First we note from 

Anderson’s result stated below ( 4.10 ) that the random rational function 

n 

Wn+l , Wi 


Rn+i{x) := 


E' 


(4.40) 


^— X — Vi 
i=l 

where (wi,..., w„+i) is distributed according to D„+i[(k)’^; a + 1], has the PDF for its zeros given 

by ji"’"^(a;,y)/ji‘’"^(y). Let us define {A|‘®(a;)}j=o,.. as specified by the recurrence ( [4.38| ) but with 
Qfo a^+ai + k, (3 q /9i throughout. Furthermore, with (wi, W2) distributed according to Bn[nk, oo + l] 
define 

Vn{x) = W2Af®(a;) + wixAft-y{x). (4.41) 


The significance of Vn{x) is seen by noting from the argument below ( 4.3C ) that with yi,...,y„ denoting 
the zeros of A^‘^{x) we have 


At{x) 


= E 


Ul 


x-yi 
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where (ui,..., is distributed according to fc], and then proceeding as in the derivation of 

( [4.31D to deduce from this the expansion 

Vn(x) Wn+l 




X 


E 


Wl 


1=1 X - y'l 


U) 


(4.42) 


where (i&i, .. .,Wn', Wn+i) is distributed according to oq + 1]. The right hand side of (4.42) is 


just the rational function (4.40), and so the PDF for its zeros, given {yi ,..., ?/„}, is 






(4.43) 


But the marginal distribution of {yi, ..., ?/„} is J^'~\y) with a = ap + ai + fc, /3 = /3i. Multiplying this 
by ( 4.43| ) shows that the joint distribution of the zeros of {Vn{x), {y)} is given by ji^''^\x,y)\a=ao- 

We note from ( ^.31 ) that if we set ai = 0, (5i = I3 q — k in the construction of {^|^®(a;)}j=o,...,n and 


then compute Vn{x) according to (4.41), the marginal distribution of the zeros of Vn{x) are given by 


ji”’‘^(x)| <»=°o . But according to (f5.3l|) and (|3.26|) the latter is identical to Jo”’‘^(x)| °=°o and so the 

01=0Q—^ '- ' ^' /3 = ^0 

marginal distribution of the zeros of Vn{x) in this case is the same as that for the zeros of A^{x). Thus 


with {wi,W 2 ) as in ( 4.41 ) we have 

A*{x) = [w 2 A*{x) +wixAf_^{x)^ 


/3o^/3Q-fc 


(4.44) 


take the limit k 


The random recurrences ( 4.38 ) and ( 4.44 ) assume definite forms if we write ap = ak, (3 q = bk and 
oo. Thus the random variables {wq\w[^\w 2 ^'^) in (4.38) crystallize to the definite 

reads 


value {{n — j + b)/d, (j — l)/d, (n — j + a)/d) where d = 2n — j —\ -\-a + b and so 
{2n-j-l + a + h)Af{x) = {n-j + a){x-l)Af_^{x) + {n-j + h)xAf_^{x) + {j-l)x{x-l)A*_ 2 {x) (4.45) 


a = a + l 
0Q>-^b~l 


(4.46) 


with initial conditions A^-j^{x) = 0, Af{x) = 1. Similarly ( 4.44 ) reads 

{a + n)A*{x) = {aA*{x)+nxAf_^{x)^ 

Using standard Jacobi polynomial recurrences we can show that the solution of ( 4.45| ) is given by 

A*{x) = (4,47) 


where the use of'indicates the monic version of the corresponding polynomial. Furthermore, (4.47) 
satisfies ( 4.46 ). The fact from ( 4.47 ) that Af[(x) = pjf' ^^( x) can be anticipated. Thus in general 

the PDF for the marginal distribution of the zeros of A^{x) is given by ji"’"^(x), and if we put ap = ak, 
Po = bk and take the limit fc ^ oo it is a known result |3^ that this PDF crystallizes at the zeros of the 
Jacobi polynomial Pi°'~^’'^~^\x). 

In Appendix D we make use of our ability to sample from Jp"’“^ to give a Monte Carlo evaluation of 
a multidimensional integral formula 0 for the bulk two-point correlation function of matrix ensembles 
with symplectic symmetry. 


5 Random matrix realizations of the Laguerre interpolating en¬ 
sembles 


Consider the Dirichlet distribution (4. Let Sn = L and scale wi,W 2 , - ■■ ,Wn-i so that Wi WijL, 
{i = l,...,n — 1). Then in the limit L —>■ oo (4.10) reduces to the product of independent gamma 
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distributions 


-(n-l) 


n—1 






Wi > 0 


(5.1) 


r(si)---r(s„_i) 

with cr = 1. If we also scale ^ ai/L (i = 1,..., n — 1), A \/L and set a„ = 1 then the limiting form 

of the random rational function ( [4.7D reads 

n—1 

i?L(A):=l + ^^^ (5.2) 


i=l 


where the Wi are distributed according to (|5.l| ) (the superscipt “L” denotes Laguerre). The distribution 
of the roots of (p]2) is given by the appropriate limiting form of the Anderson density (|4.3|), in accordance 


with Anderson’s result stated below (4.10). Let us take note of the explicit form. 


Corollary 3. Consider the random rational function {5.i) with the coefficients Wi,, Wn-i distributed 


according to (5.1). This has exactly n—1 roots, which since the Wi are positive, are real. For given {oi} 
these roots have the PDF 


1 


r(si) • • •r(s„_i) 


=> H 7 = 1 (^3 


n 


(A.-A,) 


(oi - 


:t n I* 


i-aj\ 


where 


l<2<j<n 


Ai > ai > A2 > Ct2 > ■ • • > An > Gn- 


1,3 = 1 


(5.3) 


(5.4) 


We remark that this PDF is implicit in the work of Evans [Q, who was studying Laguerre analogues of 
the Selberg integral using the method of Anderson. Also, special cases of Corollary |3| are known from [|^. 

A matrix structure for which the eigenvalue condition reduces to the calculation of the roots of 
(^.2[) (with n replaced by n + 1 for convenience) is easy to specify. One approach would be to consider 


the appropriate limiting form of (4.4). Alternatively we can write down the required matrix and check 
directly that it has the sought property. Thus let A be a real symmetric (complex Hermitian) matrix 
with eigenvalues ai > 02 > • • • > a„ of multiplicities mi,... respectively. Let Ai be a vector of 
independent real standard Gaussians (complex Gaussians) having the same number of rows as A, and 
consider the matrix 

M = A + bxx\ b > 0. (5.5) 


A simple calculation along the lines of Lemma and the derivation of (4.4) shows that M has eigenvalues 
Gi with multiplicities mi — 1, and the remaining n eigenvalues are given by the zeros of the rational 


function (^) (with n n + 1). The wt in the latter are distributed according to (^), with a = 2b/(3 
and and Si = j3mi/2 where /3 = 1 (real case), (3 = 2 (complex case). Thus (5.3) with at 1 -^ (3ai/2b, 
Xi I—> (3\i/2b, gives the eigenvalue PDF of M. 

The eigenvalue PDF of (|5.5[) in the complex case with all eigenvalues distinct can be derived in a 


different way which has the advantage of applying to a more general matrix structure. This is given in 
Appendix E. 


5.1 Construction of interpolating Laguerre matrix ensembles 

We can use knowledge of the eigenvalue PDF of (^^) to construct matrix ensembles which have as 
their eigenvalue PDF ( |1.l| ) and ( [1.21 ). For the parameter dependent PDF ( |1 . ![ ) we begin by recalling 
(see e.g. |^) that 0 is realized as the eigenvalue distribution of the matrix ensemble of 4n x 4n 
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antisymmetric Gaussian random matrices, in which the elements are pure imaginary numbers with each 
2x2 block having a real quaternion structure, and one changes variables Af i—> A^. Such matrices are 
equivalent to block matrices of the form 


02n x2n X2n X 2n 
^2nx2n ^2nx2n 


(5.6) 


where X is an antisymmetric Gaussian complex matrix, real and imaginary parts having variance 1/2, 
so X^X has eigenvalue PDF (1.8) with each eigenvalue doubly degenerate. 


Theorem 3. Let X be a 2n x 2n antisymmetrie complex Gaussian matrix, and let x be a 2nx 1 complex 
Gaussian vector, where the real and imaginary parts of the complex Gaussians have varianee 1/2. The 
random matrix 

M = X'^X + (5.7) 

has eigenvalue PDF 0 with A = 1 — 2/6. 


Proof. The eigenvalues of have multiplicity 2, and thus according to the result noted below (5.5) 
the eigenvalues of M consist of the eigenvalues of X^X with multiplicity 1 (j/i,..., say), as well as n 
additional eigenvalues xi,..., Xn which must satisfy the interlacing condition 


Xi > yi > X2 > y2 > ■ ■ ■ > Xn > yn 


(5.8) 


Setting = 2, /3 = 2 in the result noted below (5.5), it follows from (5.3) that the conditional distribution 
of the x’s given the y’s is proportional to 


c-Ur=iGi-yi)/f> 


n 


(Xi - Xj) 


n 


l<i<j<n ij = l 




Multiplying this by (1.8) (with the a;’s relabelled y's) gives that the eigenvalue PDF of M is proportional 
to 

n 

g-ElhiGi + Gi-yil/b) - y^) Y[ \Xi - yj\, 

and relabelling the eigenvalues gives the desired result. □ 

By construction 

even(M) = LSE„|a=o 


where even(M) refers to the distribution of the even labelled eigenvalues of the random matrix M. 
The result holds independent of the parameter 6 in (|5.7| ). Because Theorem |3| also tells us that when 
6 = 2 the random matrix M has the same PDF as LOE 2 „|a=o matrices, we thus have a matrix theoretic 
understanding of the relation 

even(LOE 2 „|a=o) = LSE„|a=o- (5.9) 

We remark that with X denoting the (2n+ 1) x 2n random matrix which is obtained from X in ( p.7D 
by adjoining an extra row '/bx, we have 


X^X = X^X + 6xf^ (5.10) 

Thus we can interpret Theorem || as applying to the square of the singular values of X. This latter 
viewpoint indicates a special property of the case 6 = 1. Then X is identical to the (2n + 1) x (2n + 1) 
version of X with the last column removed. But it is a standard result that the singular values of a matrix 
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interleave with those of a matrix obtained by removing a single row or column (starting with the largest 
singular value of the larger matrix). Now we know that a (2n + 1) x (2n + 1) antisymmetric complex 
Gaussian matrix is such that X^X has one zero eigenvalue, and n doubly degenerate eigenvalues, the 
latter having PDF LSE„|a= 2 - The interlacing property now implies 


odd(M|b=i) =LSE„|a= 2 . (5.11) 

As done in [^, this can be checked directly from the PDF dO)- 

Consider next the PDF ( |l.2| ). We know that as A ^ — oo this reduces to ( |l. 1 1[ ) , which is the matrix 
ensemble LUE„|a=o- It is well known (see e.g. [^) that the eigenvalue PDF for the LUEn, is realized 
by matrices of the form X"^X where A1 is an n x n matrix with i.i.d. complex Gaussian entries. Let us 
replace each complex element a; + zy of X by its 2 x 2 real matrix representation ( l.lg ) and denote the 
corresponding 2n x 2n matrix by X. The eigenvalues of X are the eigenvalues of X except that in X 
each has multiplicity 2. 


Theorem 4. Let X be the 2n x 2n real matrix constructed from the n x n complex Gaussian matrix as 
specified above, and let x denote a 2nx 1 real Gaussian vector, where in the Gaussians each independent 
part has variance 1/2. The Gaussian random matrix 


M = X^X + bxxF 


has eigenvalue PDF {l.i) with A = 1 — 2/6. 

Proof. Arguing as in the proof of Theorem || we see that the eigenvalues of M consist of the eigenvalues 
of X^X with multiplicity 1, yi,... ,yn say, as well as n additional eigenvalues xi,... ,Xn which must 
satisfy the interlacing condition (^). Making use of the result noted below (^) with mz = 2, /3 = 1, it 
follows from (^) that the conditional distribution of the cc’s given the y’s is proportional to 


n 


ivi -Vj) 


l<2<j<n 

Relabelling the x’s in dU) by y's, and forming the product with the conditional distribution shows that 
the eigenvalue PDF of M is proportional to 


^-T.l:=i(.Vi + {xi-yi)/b) 




□ 


l<i<_7<n 

and relabelling this coincides with o with A = 1 — 2/6 therein. 

Analogous remarks made after Theorem ^ also apply to Theorem ^ Thus by construction 

even(M) = LUE„|a=o, 

and because when b — 2 the random matrix M has the same eigenvalue PDF as LOE„|Q=oULOE„|a=o 
matrices, we thus have a matrix theoretic understanding of the relation 

even(LOE„|a=o U LOE„|a=o) = LUE„|a=o- 


Also, the equation (5.10) holds with X replacing X. Here the matrix X in the case 6 = 1 is equivalent to 
the (2n + 2) x 2n version of X with the last row removed. Now before removing the row, such Gaussian 
matrices multiplied by their transpose have a doubly degenerate zero eigenvalue and n doubly degenerate 
eigenvalues with PDF LUE„|a=i. Arguing as in the derivation of (^.ll|) we therefore conclude 


odd(M|b=i) =LUE„U=i, 


(5.12) 


which like (5.11) can be verified by direct integration of the PDF (|l.2|) with A = — 1 
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5.2 Laguerre limit of the three term recurrences 


The joint PDF ( ^.24 ) has a well defined Laguerre limit, specified by changing variables Xi ^ XijL, 
Vi ^ Ui/L, setting (3 = L/b, /3i = L/bi and taking the limit L ^ oo. This gives 


- (n,n —1) 


{x,y) := 


_ i _ t _TT 

r(i + a)(r(fc))"-i 

i{a + ai + /c, fc; bbi/{b + bi)) 


.a-Xilb 


n 


l< 2 <j<n 






n \y3-y^'W'W\x3-Vi 

2 = 1 j—1 


k-1 


(5.13) 


where 


Wn{a, k;b) = / dxi 
Jo 


dXn 


a -xilb 


X] e 


\Xk-' 


_ ^n(o+l+(n—l)fc) 

J = 1 


1=1 l<j<k<n 

T{l + kj)r{l + a + kij-l)) 


T{l + k) 

I — 

and the a;’s and y’s are interlaced according to 

Xi > 2/1 > X2 > 2/2 > ■ ■ • > 2/n-l > Xn>0. 

Also of interest is the Laguerre limit of the marginal distribution 

1 


(5.14) 


f dyi ■ ■ ■ dyn-i ^Kx,y) 

JR 


n \xk-x,\>^=:Li-’-Hx) (5.15) 

l/bi=o Wn{a,k,b) l<j<k<n 

where R refers to the region ( ^.14| ). 

Sampling from Lo"’" (x, y) and {x) can be undertaken by taking the Laguerre limit of the three 

term recurrences in SectionConsider first the recurrence (4.38) determining the polynomial with zeros 
realizing the PDF J^’'\x). The Laguerre limit is obtained by scaling x ^ x/L, w[^'^ v[^'^ jL, 1 X 2 '^ 

V 2 ^ jL, = 1, where the v[^\v 2 '^ are distributed according to the gamma distributions r[6; {j — l)fc], 
r[&; (n — j)k + ao + 1 ] respectively (here the notation r[cr; s] refers to the density function proportional 
to With v^^\v 2 '^ so specified, and introducing the further scaling Af{x) = L~^Bf{x) we 

see that the Laguerre limit of (4.38|) reads 


Bf{x) = {x - V 2 ’)Bf_^{x) - xvY’Bf_ 2 {x). 


.i3)\ TJ# 


U) D# 


(5.16) 


This recurrence is to be solved subject to the initial conditions B'^^{x) = 0 and Bf (x) = 1. The zeros of 
Bfj;{x) are then distributed according to Lo"’'^(a;)|a=ao- 

To sample from Lo"’" ^^(a:, y) we take an appropriate Laguerre limit of the procedure to sample from 
Jo”’"~i)(a:, 2 /) detailed below (t4.38|) . Thus we use ( 5.16| ) to first compute {i3j^(a;)}j=o,...,rt-i where B^{x) 
refers to B'f {x) with parameters ao ao + cti, 1/6 i—> 1/6 + I/ 61 . We then form the random polynomial 


S„(x) = {x- vi^'>)B*_j^ix) - xv'f^B*_ 2 {x) (5.17) 

with distributed according to (r[l/6; (n — l)fc], r[l/6; ao + 1]). The zeros of (i?„(a;), Bl^_^{y)) 

then have the joint PDF Lo"’” ^\x,y)- 

A recurrence to sample from L^’~\x)\ a=a-k(,n-i)-i has been given by Dumitriu and Edelman |]l^ , 
as a corollary of their construction of a random tridiagonal matrix with this eigenvalue PDF. Denote 
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by Xp the gamma distribution r[2;p/2], and let Xp denote the square root of a random variable with 
distribution Xp- It was shown in that the symmetric n x n random tridiagonal matrix 

G’n 1 

1 ^n—1 ^n—2 

bn-2 an-2 &n-3 


02 

bi 


bi 

ai 


(5.18) 


with the distribution of the elements specified by 

On ^ X2ai ^ Xki d" X2a-fe(n-i)> ^ XkiX2a-k{n-i-l) (* = "O — 1, . . . , 1) 
has eigenvalue PDF given by a=a-k{n-i)-i. The random recurrence now follows from the fact 


that in general the characteristic polynomial of the bottom right k x k submatrix of (5.18) satisfies the 
three term recurrence 

Pj{x) = {x- aj)Pj-i{x) - b]_^Pj- 2 {x), (5.19) 

subject to the initial conditions P-i{x) = 0, Po{x) = 1. Note that ( 5.19| ) differs from ( 5.16| ). 

The distribution Je^’^\x,y) as given by (|3.29| ) also has a well defined Laguerre limit, obtained by 
writing Xi (1 - Xn+i-i/L), yi {1 - yn+i-t/L), setting a = L/b, ai = L/bi, f3i = ai and taking 
L ^ oo. This gives 

1 1 


L("’")(x,y) = 


r(l + Q!)(r(A:))"- Wniai^k, (6 + bi)/hhi) 

n 7 

11 \^3 - Xi\\yj - yi\ \ 

i—1 

where the x’s and y's are interlaced according to 

xi > yi > X 2 > 2/2 > • • • > a;„ > > 0. 

We make note too of the Laguerre limit of ( |3.31 ), 

1 


(5.20) 


(5.21) 


dyi ■ ■ ■dynL^J"’'^\x,y) 


IR 


^ai+k^-Xilb 

l/bi=0 Wn{ai + k, k, b) 




(5.22) 

where R denotes the region ( ^.21 ). 

To sample from L^’'^\x,y), we note that the Laguerre limit of {A|^'^(x)}j=o,...,n, used in ( 4.41 ) to 
sample from (x, ?/), is given by {B]^'*(x)}j=o,...,ra where the B^‘^{x) are specified by the recurrence 

(^.16) but with 1/6 1 -^ 1/6+ l/6i, ao i—> ai in the specification of the distribution of (up\ t> 2 '^^). Then 


taking the Laguerre limit of (4.41) tells us if we define 

Un{x) =U 2 B*^{x) +UiBfZi{x), UiPU2 = l (5.23) 

where ui is distributed according to r[6;nfc], then the joint distribution of the zeros {f7„(x), i3^®(y)} 
is given by Li^’^\x,y). We remark that according to ( p.22 ) the marginal distribution of Unix) in the 
case l/6i = 0 is given by L^’'\x,y). But i?)f®(x) with l/6i = 0, ao ai is the same as B^(x) with 
ao 1 -^ ai and thus has the PDF for its zeros given by (^.15|) with a = ai, and this in turn is identical to 


Li”’‘^(x) 


Thus we have the Laguerre limit of (4.44), 

B*{x) = (u 2 B*{x) + UixB*_-^ix)^ 


ai I—*^Q;i —k 
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6 Gaussian interpolating ensembles 

The random rational functions o and ( |5.2[ ) have a counterpart 

n 

R^{X) := A + 




a-i — y 


( 6 . 1 ) 


where Wj is distributed according to r[l;sj], which is closely related to the Gaussian ensembles. Thus 
implicit in the work of Evans on Gaussian analogues of the Selberg integral according to the method 
of Anderson is the following result for the PDF of the zeros of R^{X). 

Proposition 8. Consider the random rational function Q. This has n+1 zeros Ai,..., A„+i restricted 
by the interlacing condition 

Ai > ai > A 2 > 02 > • • • > On > Ara+l 


and the further requirement that 


n+1 


1^1 1^1 

Subject to (^) and B the PDF of the zeros of B is given by 


( 6 . 2 ) 

(6.3) 


1 Y\l<j<k<n+li^3 Afc) 


n+1 


r(si) • • • r(s„) Ok) 






- Ob 


1=1 P=1 

Also of interest is the random rational function 




1=1 


1=1 


R^(X) := A - Wo + , , 

•“ ax — A 


(6.5) 


1=1 


where the Wj (j = 1,... ,n) are distributed as in (6.1) while wq is distributed according to N[0,1]. The 
PDF for the zeros of ( |^ ) is readily deduced from Proposition |[ 

Corollary 4. The zeros of the random rational function (6.5) have PDF (B multiplied by XjyfTK, 
except that the condition B is no longer required. 

Proof. The rational function i?’^(A) results from R(X) by making the replacements A 1 -^ A — wq, 


Oj — Wq. We then have that 


^ n+1 ffc /t-pj. ffc 


n+1 


exp 


e^o/2 


exp 


n+1 


n+1 


1=1 


1=1 


j=i 


1=1 


1=1 


1=1 


1=1 1=1 

Multiplying this by (the distribution of wq) and integrating over wq eliminates the delta 

function (and thus the restriction ( |6.3| )) but leaves all other terms unchanged. □ 

We remark that the random rational function (|6.5|) can be derived as a limit of the random rational 


1 

IF 


function ( 4.21 ) with n n+1 and (wq, ..., Wn\ w„+i) distributed according to D„+ 2 [q!/ 2 , si,..., Sn] a/2]. 
Thus if we write a = Lf and take L —> 00 then the marginal distribution of wq and w„+i have the 
asymptotic form A + ^N[0,1] while the Wi (i = 1 ,... ,n) have to leading order the marginal distribution 
r[l; Si]. It then follows from ( |4.21 ) with x +1 — ^) and yi 1 —s- i(l — ^) that we have 

j-Rn+ 2 {x) ^ ++A). (6.6) 

Z L—*oo 
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The random rational functions 0) and (|6.5|) occur in two closely related eigenvalue problems (see 
e.g. 1^). Thus let A be a real symmetric (complex Hermitian) matrix with eigenvalues oi > 02 > • • • > a„ 
of multiplicities From A form a random matrix M of one extra column and one extra 

row by bordering A by a constant y/b times a vector of independent real standard Gaussians (complex 
Gaussians) as the final column, and the Hermitian conjugate of this as the hnal row (therefore in both 
the real (/3 = 1) and complex case (/3 = 2) we require the final entry of the vector to be real; let it have 
distribution N[0,-\/2//3]). Thus if the final column of A is number n*, then 

[M]i^j=A, [M],^„-+i = = Vb[x]i, {l<i,j<n*) ~ iV[0, (6-7) 


where here the symbol ^ denotes ‘has distribution’. A straight forward calculation shows that M has 
eigenvalues ai with multiplicities — 1, and n + 1 further eigenvalues given by the zeros of the rational 
function (^) with wq (wj) distributed according to N[0, (T[2b/(3, firrij/2]). It follows by scaling 

(H) that if we choose c = 2b then the eigenvalue PDF of M is given by ( |6.4| ) with Xj f3/2bXj, 

If in the prescription (HI) we choose [M]n*+i^n*+i = 0 we find that M has eigenvalues Oi with 
multiplicities rrii — 1, and n + I further eigenvalues given by the zeros of the random rational function 
( |6.l[ ) with Wi as specified in the above paragraph. Note that the condition (^(^) then has the interpretation 
as the statement that Tr(A) = Tr(M). 


6.1 Construction of Gaussian interpolating matrix ensembles 

Following a strategy analogous to that used in the construction of random matrices realizing the Jacobi 
and Laguerre interpolating ensembles, we can use the eigenvalue problem relating to (^) to construct 
random matrices with eigenvalue PDFs realizing certain Gaussian interpolating ensembles. In particular 
we can construct random matrices with eigenvalue PDF of the form 


n+1 

n- 

2=1 


-cix\l2 


n ixi-Xj)Y[e iVi-yj), 




2=1 


l<2<j<n 


where 


xi > yi > ■ ■ ■ > yn > Xn+l, 


which with C 2 = 0 reduces to (2.48), and the eigenvalue PDF 


( 6 . 8 ) 


(6.9) 


n+1 




3 cia;2j_i/2 


n- 


-C2xlj2 


2=1 


n 

l< 2 <j'< 2 n+l 


where 


Xi > X2 > ■ ■ ■ > X2n+1- 


( 6 . 10 ) 


( 6 . 11 ) 


To obtain (6.8) we choose the matrix A in (6.7) to be an n x n member of the GUE (see e.g. [Q 
for the precise definition of such matrices), and we extend A so specified to a 2n x 2n real matrix by 
replacing each complex element by its 2 x 2 real matrix representation ( |4.18| ). Following the strategies 
of the proofs of Theorems ^ and ^ it follows that M as specified by (H) with n* = 2n, /? = 1 has the 
n eigenvalues of A with multiplicity 1, yi,..., say, and a further n + 1 eigenvalues xi,..., Xn+i say 
interlaced according to (^^). Furthermore it follows that the joint eigenvalue PDF is given by (^) with 


2 &’ 




( 6 . 12 ) 
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The above construction gives 

even(M) = GUE„. (6.13) 

Nowwith b = 1/2 we have ci = C 2 = 1 and we recognize ( ^ ) as the eigenvalue PDF for GOE„+iUGOE„ 
(see e.g. Q). Thus we have a matrix theoretic understanding of the identity fl^ 

even(GOE„+i U GOE„) = GUE„. (6.14) 


We note too that with & = 1/4 the matrix M coincides with the upper left (2n + 1) x (2n + 1) block 
of the real matrix representation of a (n + 1) x (n + 1) GUE matrix. By an argument analogous to the 


derivation of (5.11) we must therefore have 


odd(M|5^i/4) — GUE„_|_i. 


(6.15) 


As 6 = 1/4 corresponds to C 2 = 0, this identity is relevant to (|2.48| ). 

For the realization of ( 6.10| ) we choose the matrix A in ( |6.7D to be a n x n member of the GSE (by 
definition — see e.g. |l^ — the elements of such matrices are real quaternions, so as a complex matrix 
A is 2n X 2n). The eigenvalues are doubly degenerate, with the independent eigenvalues j/i, ... say 
having distribution 

n 

n (yj-ykf- 

1—1 

Here we follow the strategy of the proofs of Theorems ^ and ^ to conclude that with this choice 
of A and n* = 2n, [3 = 2 has the n eigenvalues of A and with multiplicity 1, a; 2 , a; 4 ,..., X 2 n say, and a 
further n + 1 eigenvalues xi,a; 3 ,... ,X 2 n+i say, interlaced according to ( 6.1l| ) and with eigenvalue PDF 
( |6.10 ) with 


1 1 
Cl = T, C 2 = — - 


2 . 


5’ b 

Since by construction 

even(M) = GSE„, 

and with b = I and thus ci = C 2 = 1 the PDF ( |6.10 ) reduces to the PDF for GOE 2 n+i, we thus have a 
matrix theoretic understanding of the relation fl^ 


even(GOE 2 „+i) = GSE„ 


Furthermore, with & = 1/2 the matrix M coincides with the upper left (2n + 1) x (2n + 1) block of the 
complex representation of a (n + 1) x (n + 1) GSE matrix, and so we must have 

odd(M|5^i/2) = GSE„_|_i. 


6.2 Gaussian limit of the three term recurrences 


The PDFs ( |6.^ ) and ( |6.10D are special cases of a limiting form of the joint PDF (3.24). Thus in ( ^.24 ) let 
us change variables Xi — ^), (| — ^), set a = [3 = aL'^, ai = /3i = aiL^ and take L oo. 

We then obtain the joint PDF 


y) •= 

{.x,y) . (r(fc))"-i 


1 


(r(fc))”-i M„_i(fc; 2 (a + ai)) 




X I I e 


n 

l<z<j<n 


Xj-x,\W_e n \y3-y^\W_W_\xo-y3 ^ 


(6.16) 


2=1 


l<i<j<n—l 


2=1 i=i 


35 
























where 


/ oo poo 

dxi--- da;„ JJ 

-OO ^ —OO t ^ ■ ^ ■ 


\2k 


l<i<j<n 


_ ^—nj2 — kn{n—l)j2 


(27r)"/2 n 
j=o 


r(i + {j + i)k) 


T{l + k) 

j — 

and the a;’s and y’s are interlaced according to 

oo > xi > yi > X 2 > y 2 > ■ ■ ■ > Un-i > Xn> -oo. 
We make note of the special marginal distribution 


f dj/i • •-dyn-i ^Hx,y) 

Jr 


1 


ai=o M„(A:;2a)n 


— 2ax^ 


n 1 ^* “ ^3 

l< 2 <j<n 


2k 


(6.17) 


(6.18) 


which is equivalent to the integration formulas (3.26) and (5.15). 

To sample from G^’~\x) we can take the Gaussian limit of the three term recurrence (4.38). First we 
note that with {wq\wi\w 2 ^) distributed as specified below ( 4.38 ), setting ao = Po = alP and taking 


L —> oo, the marginal distributions of and W 2 have the asymptotic form i 

J. 

I2a 

'"j 


:N[0, 


■ 7=1 while ic, 

/la -* 


(j) 


has the leading order marginal distribution {j — l)fc] (c.f. the statements above ( | 6 . 6 | )). Thus by 

also writing x ^(1 ~ f)’ (^) {—‘JL)~3Cf (a:), we see that in the Gaussian limit ( 4.38|) reduces to 


Cf{x) = (x - r)Cf_i{x) - ^'‘Cf_ 2 {x) 


(6.19) 


where r has distribution N[0, while s^3 1 ) has distribution r[^; (j — l)fc]. With the initial conditions 
G)f;^(a:) = 0, Cf{x) = 1, we have that the zeros of Cf{x) have PDF G^"\x). The recurrence (6.19) has 
the structure (5.19) and thus can be viewed as specifying the characteristic polynomial for a corresponding 
random tridiagonal matrix ( |).l^ ). In fact this is precisely the random tridiagonal matrix found by 
Dumitriu and Edelman fl^ and shown to have eigenvalue PDF given by G^’~\x). 

To sample from Go"’" ^\x,y), the Gaussian limit of the procedure to sample from Jo"’" ^\x,y) 
detailed below ( 4.38| ). Thus we use ( |6.19| ) to generate {Gj^(a:)}j=o,....n-i where Cf{x) refers to Gf{x) 
with parameter a 1 —> a + oi. We then form the random polynomial 

Cn{x) = {x- r)cf_^{x) - s^^-^'>C*_2{x) 


with (r,s^" ^)) distributed according to (N[0,F)^; (n — l)fc]). The PDF G^"’" ^\x,y) is then 
realized by the zeros of (G„(a;), C^_i{y)). Equivalently we can realize the PDF in terms of the eigenvalues 
of a random n x n tridiagonal matrix and its lower right (n — 1 ) x (n — 1 ) submatrix. 

Theorem 5. Consider the symmetric tridiagonal matrix Let the elements be random with dis¬ 

tributions 


a* ~ N[0, 




FI 


1 


2(a + fli) 


\/2(a + oi) 

N[0,^], bl_,r^T[^;{n-l)k]. 


;{i-l)k] (i = l,...,n- 1 ) 




'2a 


The joint distribution of the eigenvalues xi,... ,Xn of this matrix, and the eigenvalues j/i,..., yn-i of the 
(n — 1) X (n — 1) bottom right submatrix, is given by Go"’" ^\x,y). 
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Appendix A 


Geometrical RSK 


The Robinson-Schensted-Knuth correspondence between non-negative integer matrices and pairs of semi¬ 
standard tableaux of the same shape has a geometrical representation relating the matrix to paths defining 
the interface of a sequence of growth models |2^ (for closely related geometrical representations see @ 
and 1^). Here will will show that coordinates specifying the paths can be used to deduce the equation 
( ^.10 ) and also the joint probability (2.4). 

First we revise the construction of 1^ , wherein each distinct n x n non-negative square matrix 
X = (rows counted from the bottom), with entries Xij weighted (1 — aibj){aibjY''\ is 

put into a one-to-one correspondence with a set of at most n non-intersecting weighted lattice paths, 
starting at (x, y) = (—(2n — 1/2), Z — 1) and finishing at (x, y) = ((2n — 1/2), Z — 1) (Z = 1, 2,...). The 
Zth member of the set — the path starting and finishing along y = — {I — \) — will be referred to as the 
level-Z path. Each level-Z path can be regarded as a pair of paths because the weights and the allowed 
steps are different depending on x < 0 or x > 0. Thus the first (second) member of the pair starts at 
X = — (2n — 1/2) (x = (2n — 1/2)) and goes either right (left) in steps of two units, or up in integer 
amounts at x = —(2n -I- 3/2 — 2j) (x = (2n -|- 3/2 — 2j) with each unit regarded as a step weighted by bj 
(aj), until it reaches x = —1/2 (x = 1/2) where both paths must have the same final y-coordinate (see 
Figure ||). 
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Figure 1: The RSK mapping from a weighted integer matrix to a set of weighted non-intersecting 
paths. The mapping is invertible and so is a bijection. 

Let yi denote the maximum height of the level-Z path, which is the displacement of this path at 
X = ±1/2. Because the paths cannot intersect we must have p-i > /22 > • • • > so y = (yi, / 22 ,..., /r„) 
forms a partition. It is a standard result (see e.g. |Q) that the total weight of all non-intersecting paths 
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of the specified type for x < 0, initially equally spaced at y = 0,..., — (n — 1) along x = — (2n — 1/2) and 
finishing at y = /ii, /j ,2 — 1,..., /J.„ — (n — 1) along x = —1/2 is given by the Schur polynomial s^(6i,..., 6„). 
Similarly the total weight of all non-intersecting paths of the specified type for x > 0, initially equally 
spaced at y = 0,..., — (n — 1) along x = {2n — 1/2) and finishing at y = fj,i, ^2 — , ^J‘n — {n — 1) along 

X = 1/2 is given by the Schur polynomial s^(ai,... ,an)- Furthermore it is another standard result (see 
e.g. Q) that each set of non-intersecting paths from x = — (2n — 1/2) to x = —1/2 is equivalent to a 
semi-standard tableau of shape y, content n, as each set of paths from x = (2n — 1/2) to x = 1/2. Thus 
for the non-intersecting lattice paths of Figure]^ we obtain the pair of tableaux 


1 

1 

2 

2 

2 

3 

2 

3 



1 

1 

1 

1 

3 

3 

2 

3 



the first corresponding to the paths x > 0 and the second the paths x < 0. In the former (latter) each 
■j j in the level-/ path is recorded as a box labelled j in the /th row of the tableau. 


occurence of aj (h 


Thus if we accept the mapping, the probability (^) is immediate (an unequal number of a and b 
weights can be achieved by simply setting some of them equal to zero and using the stability property of 
the Schur polynomial, s^(ai,..., a„_i, 0) = s^(ai,... ,a„_i)). 

To derive (|]J) and (|2.1Cl|) we must investigate the details of the mapping, which takes the form 
of a cascade of polynuclear growth models (see Figure |^. We rotate the non-negative matrix X = 
45° anti-clockwise and label the horizontal rows of the rotated matrix by t = 1, 2 ,..., 2 n— 1 
and the vertical columns byx = 0,±l,...,±(n—1) where x = 0 corresponds to the diagonal i = j oi X 
(recall that the rows are being counted from the bottom). The entries Xij in the matrix for successive t 
values (t = i+ j — 1) are heights of weighted ‘nucleation events’ — columns of unit width and height Xij 
centred about the corresponding x-coordinate which are placed on top of the profile formed by earlier 
nucleation events and their growth and weighted by (0^6^)“••■» (in addition the matrix has a normalization 
weighting of aibj) independent of the entries). Thus at t = 1 there is a nucleation event at 

X = 0 which consists of a column of width 1, height xn and weight (0161)’”'^'^ marked on the line at y = 0 
in the xy-plane. In general, as 11—> t-l- 1 the profile of all nucleation events so far recorded is to ‘grow’ one 
unit in the —x direction and one unit in the +x direction. Thus in going from t = 1 to i = 2 the nucleation 
event centred at x = 0 of height xn now has width 3 units. On top of this profile, centred at x = —1 
and X = 1 nucleation events of unit width and height X2 i,xi 2 and weight {a 2 bi)^^^, (0162)’”^^ respectively 
are then drawn. In now going from t = 2 to t = 3 this new profile is to grow one unit to the left and 
one unit to the right. In so doing we see that an overlap of width one unit and height min(x2 i,xi 2) 
will occur. This overlap is ignored in the first diagram (profile on y = 0), and recorded instead as a 
profile on the line immediately below (here y = —1). The process is repeated with these rules until the 
nucleation event of height x„„, weight (a„/)„)‘’’"" at i = 2n — 1 has been recorded above x = 0 on the 
first diagram. In this way we obtain the sought mapping from a weighted non-negative integer matrix to 
weighted non-intersecting lattice paths. The mapping is easily seen to be invertible, and so is a bijection. 
With each path considered as a pair of paths depending on whether x < 0 or x > 0, and then the set of 
paths for x < 0 and the set of paths for x > 0 recorded as a pair of semi-standard tableaux, this gives 
the same correspondence between non-negative integer matrices and pairs of semi-standard tableaux of 
the same shape as the Robinson-Schensted-Knuth algorithm (this last point follows because, as noted in 
the above algorithm can be viewed as a graphical presentation of the matrix-ball construction of 
Fulton 11^ , which has been shown to give the RSK correspondence). 

Of crucial interest to us is the sequence of maximum displacements A/(711,712) of the level-/ path 
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obtained by applying the growth process to the truncation say of the matrix X to the first ni rows 

and 712 columns, extended to a square matrix of dimension max(ni, 772) x max(ni, 712) by appending rows 
of zeros to the top (or columns of zeros to the right, as appropriate). Consider first the level-1 path. It 
follows from the rules of the growth process that Ai( 771 , 772 ) results by adding Xn-i,n 2 to the maximum of 
the height at a: = —1 and the height at a: = 1 in the previous time step {hi{ni, 712 — 1) and hi(ni — 1, 772) 
respectively say; see the change in the height at the origin in going from t = 4 to t = 5 in Figure Thus 

Ai ( 771 , 772 ) = max (771,772 - l),hi{ni - 1 , 772 )) -|- Xni,n2- 

But again from the rules of the growth process 

/ll( 77 l, 772 - 1) = Al ( 771,772 - 1), /ll (tIi - 1, 772) = Al ( 77 l - 1 , 772) (A.l) 


since the nucleation events in column 772, and rows 1, 2, ..., 771 — 1 from the bottom cannot contribute to 
hi { ni , 772 ~ 1), and similarly the nucleation events in row 771 from the bottom and columns 1, 2, ..., 772 — 1 
cannot contribute to hi{ni — 1,772) (see Figure Therefore we obtain the recurrence 

Ai( 77 i, 772 ) = max (Al (771, 772 - 1 ),Ai(? 7 i - 1,772)) +Xnj_,n2 (A. 2 ) 


which with the boundary condition 


Ai(0,j) = Ai(7,0) = 0 


uniquely specifies {Xi(i, One observes that L{ni,n 2 ) as specified by ( 2 . 2 ) satisfies the very 

same recurrence e, and indeed one of the primary relations in the RSK correspondence is 


Ai(ni,n2) = L(ni,n2). 


(A.3) 


0 

0 2 
0 0 0 
2 1 
1 


0 


0 

0 0 
1 0 0 
2 1 


Figure 2: The path diagram for the matrix of Figure with the final row (column) set equal to 
zero. This can be deduced from the t = 4 diagram of Figure Q with 03 ( 63 ) set equal to zero. 
The maximum heights therefore coincide with the maximum height at x = 1/2 {x = —1/2) in 
this diagram. 


For the maximum displacements Xi{ni,n 2 ) of the level-Z path, I > 1, we see from the derivation of 
Ai( 77 l, 772 ) = max ^Ai (771, 772 - l),Ai( 77 i - 1 , 772 )) + X^n^.nl 


(A.2) that 


where Xn^ is the height of an overlap event (if any) which occurs in the growth of the nucleation events 
corresponding to x„_i_„ 2 _i or Xni-i,n 2 - We remark that [xfj defines the Ith member of the 

sequence of matrices in Fulton’s matrix ball construction . The rules of the growth process give 


,n2 


min (|/7i_ 1(771, 772 - l),/7i_l(77l - 1 , 772 )) - hl^l{ni - 1,772 - 1) 
min ( Az_i( 77 i ,772 - 1), Ai_l(77l - 1 , 772 )) - Ai_l(77l - 1,772 - 1) 
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(see the nucleation events created in level 2 in going from t = 2 to t = 3, and going from t = 4 to t = 5 
in Figure and so for Z > 1 


A;(ni,n2) 


max 712 - 1), A;(7ii - 1 , 712 )^ 

+ min ^A;_i(7ii,n2 - 1), Ai_i(ni - 1 , 712 )^ - - 1 ,772 - 1), 


(A.4) 


which with the boundary condition 

Ai(0,j) = Ai(i,0) = 0 

and knowledge of {Xi{i, j)}ij=i^,,,^n from e uniquely specifies 


(A.5) 

For recurrences 

closely related to (|A.4| ), see |30|] . 

We have defined A; ( 711 , 772 ) as the maximum height of the level-Z path resulting from applying the 
growth process to the truncation „2 of the original matrix X. But there is another equally important 
interpretation of A;(77i, 772 ). Thus consider ni, 77,2 as fixed and apply the growth process to Xni,n 2 - Then 


we can see, arguing as in the justification of the equalities (A.l), that Xi{ni,j) is equal to the displacement 
of the level-Z path at x = —2n* — 3/2 + 2j, while A; (7,77*) is equal to the displacement of the level-Z path 
at X = 2n* + 3/2 — 2i, where n* = max(r7i, 772 ). It follows immediately from this interpretation that (see 
Figure ||) 


and furthermore 


= 0 for Z > max(77i,772) 

XiiiJ) > Xiii,j - 1) > Ai+i(7, j) 
XiiiJ) > Xi{i - l,j) > Xi+i{i,j), 


^ (Ai(77,j) - Xi{n,j - 1)) 

n 

^ (^Xi{i,n) - Xi{i - 1 , 77 )) 

Z=1 


n 

i=i 


(A.6) 


(A.7) 


Ar+i(b/) 


Ai(7,j) 

Ai(7 - l,i) 


Figure 3: Graphical demonstration of the final inequality in ( |A.6| ), where here all marked dis¬ 
placements are with respect to the line y = —{I — 1), and A/^_i(i,j) := Az+i(7,j) -|- 1. 


We are now in a position to derive the joint probability (2.4). Let Xn^^n^+i he an integer matrix 
mapping to a pair of semi-standard tableaux of shape y. under the RSK correspondence. Now each part 
yi of y is equal to the maximum displacement yi of the level-Z path, and so from the above discussion 
yi = A; ( 771 ,772-1-1). Similarly, with k denoting the shape of the pair of tableaux resulting from applying the 
RSK correspondence to the truncation Xn^^n^ of Xni,n 2 +i obtained by deleting the rightmost column, the 


above discussion shows ki = A; (771,772). The relations ( A.6 ) immediately give (|2.5| ) in the case 772 > 77 i 
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and (^) for n 2 < ni. Furthermore, from the geometrical RSK mapping illustrated in Figure 1^, the 
joint probability is seen to be equal to ~ o-ibj) (the normalization weighting), times the 

total weight of all non-intersecting paths from (2n* -|- 1/2, ^ — 1) to (1/2, Z — 1) (Z = l,...,n*) with 
n* = max(ni, 712 + 1), weighting Uj for each unit up step at a; = 2n — 3/2 -1- 2j, {j = 1,..., ni), times the 
total weight of all non-intersecting paths from {—2n* — 1/2, Z — 1) to (—3/2, Z — 1), weighting bj for each 
unit up step at a; = 2n* — 3/2 -|- 2j (j = 1,... ,ni), times the step weight bn^+i raised to the power of 
the difference in the maximum height of the level-Z path for a; > 0 {^i) and the level-Z path for x < 0 (k;) 
summed over Z. Writing the total weights of the paths in terms of Schur polynomials we see that ( ^.4[ ) 
results for 712 > tii, and (2.4) modified by (H) results for 772 < n-i. 

To derive ( 2.10 ) we first note that in the special case X = is symmetric about i = j, 

we must have \i{i,j) = Hence it follows from ([A.2|) and (|A.4D that then 


= Xi{i,i - 1) + Xi^i 

Xi{i,i) = Xi{i,i - 1) + Xi-i{i,i - 1) - Xi-i{i - l,i - 1), Z > 1. 
Forming appropriate linear combinations of these equations shows 

i i—1 


^(- 1 )' ^Xi{i,i) ^Xi{i- l,i- 1) = Xi^i, 

1^1 1^1 

and summing this equation over i from 1 to n gives ([2.10). 


Continuous RSK 

In the above description of the RSK correspondence a weighted non-negative integer matrix X = 
iias been mapped bijectively to a set of non-intersecting weighted lattice paths. Because 
for given maximum displacements /ri,... ,/7„ the total weight of the non-intersecting paths given by a 
product of Schur polynomials of the same index /x, this allows the probability that the integer matrix 


matrix maps to such lattice paths to be specified by (2.3). 

Also of interest is the case when the matrix X consists of non-negative real valued random variables 
Xij distributed according to the exponential distribution 


Prixij G [y,y + dy]) = {ai + Pj)e y > 0. 


(A.8) 


Using this distribution to define a probability measure on X, we see (as noted in [|4)) that the RSK 
correspondence gives a bijective mapping to a set of non-intersecting paths with a certain probability 
measure. The description of the lattice paths differs in some details to the discrete case. First, their steps 
are continuous in the y-direction and discrete in the ^-direction. All paths start along y — 0, with the 
level-Z path starting at {x, y) = {—{2.n + ‘i/2 — 2Z), 0) and finishing at {x, y) = {2n -|- 3/2 — 21, 0). Vertical 
steps can occur at x = —{2n + 3/2 — 2j) {x = 2n + 3/2 — 2j),j = Z,..., n, with the constraints that the 
non-zero height of the level-Z path is greater than that of the level-(Z -|- 1) path, and the height of each 
level-Z path is weakly increasing going from x = —(2n — 1/2) to x = —1/2 and weakly decreasing going 
from X = 1/2 to X = 271 — 1/2. 

We are interested in the probability density that the level-Z paths will have maximum displacement 


yi, I = 1,... ,n. By the choice of the exponential distribution (A.8), the RSK correspondence shows that 
at X = —{2n + 3/2 — 2j) {x = 2n + 3/2 — 2j) the vertical increment of each level-Z path with Z < j is a 
random variable proportional to (e““^^) — the normalization n"j=i (®j +/3i) is taken as an overall 
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factor — which is conditioned so that the paths for levels 1, 2,..., j do not intersect, and furthermore at 
level-/ the sum of the vertical increments for both a:: > 0 and a; < 0 is equal to yi. 

The total weight of a single path with vertical increments of length Vj at a; = —(2n -I- 3/2 — 2j) 
(j =/,..., n), weighted by and constrained so that = Vk is given by 



poo . n n 

/ dSn ^(j/fc - X! = X! 

j=l 3=1 


o-djVk 


J—L fi^j 




The extension of the Karlin-McGregor theorem used in gives that the corresponding total weight of 
the set of continuous non-intersecting paths for a: < 0 is det[ui{{Pj}j-i^,,,^n',yk)]k,i=i,...,n- Similarly for 
a: > 0 the total weight of the set of continuous paths is det[ui{{aj}j-i^,,,^n',yk)]k,i=i,...,n- Thus the sought 
probability density is 


(cTj-t- fdj') det[zi/, y/e)]det[w/, y/j;)]..,n■ 

jj=i 


In the special case 

ai = a + {i — 1), Pj = d + (j — l)c 

we have the simplifications 

g-(a+(i-l))yfc 




-(1-e ^'=)" \ ui{{Pj}j=i,...,n;yk) = 


^-{a+{l-l)c)vk 


-( 1 - 


c'j'n-l 


and, after making use of the Vandermonde determinant identity, we obtain the corresponding simplifica¬ 
tion of (A.9) 


r(a-b g-b (j - l)c-f n) (^+5) 


n r(a-ba-b{/- l)c)c-J' ir2(/)' 


n 


-cyj 


— e 




(A.IO) 


l< 2 <j<n 


Appendix B 

The summation identity ( |3.17|) 

Substituting ( |3.13 ) in (|3.17f) shows that we have the normalization identity 




(B.l) 


where the summation over k is over the region ( 3.14 ). Making use of (^) and ( |3.7| ) we find that with 

(B.2) 


^ = g^ := Xn+l-^-.= 


(B.l) can be rewritten to read 


n— 1 


En Vj n (y^-y 3 ) n im/^^:Q)k-liqyi/xJ+l■,q)k-l 

y 3 —^ 1 

{{q-,q)k-iy 


{q]q)kn-i 


{xi - xj){qxi/xj;q)k-i{qxj/xi)k-i 


where the summation over y is over the regions 

yi = Xi, qxi, q^Xi, q^x^,x^+l {i = 1,... ,n - 1). 


(B.3) 


(B.4) 
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This can be viewed as a multidimensional q-integral. It is the special case si = • • • = s„ = fc of Evans 
[ p^ g-integral 

n—1 

En Vj iVi Vj) iQyj/^ijQ)Bi-iiQyi/xj+i;q)sj+i-i 

= TT (^Xi - Xj){qx^/xj-,q)s^-iiqXj/xi)si-i, (B.5) 

and it is also the special case = 0 of Okounkov’s g-integral formula for the Macdonald polynomials 


( n (Xi - Xj){qxi/xj;q)k-i{qxj/xi)k-i^ ^ .. ., q, t) j/j iVi - Vj) 


l< 2 <j<n 


j—1 


^ n i<iyj/xz;q)k-i{qyz/xj+i;q)k-i = P^{xi,...,Xn-,q,t). (B.6) 

i<i<j<n-i {q,q)kn-i 


This latter observation allows us to write (B.6) in a structured form from which a simple derivation 
in the general v case follows. Recalling (^) and (^), and introducing the notation from 

[ p^ pg. 331] to denote the evaluation map on polynomials in n-variables ((n — l)-variables) which sets 
Xi = (j/j = we see that (|B.6D can be rewritten 


(B.7) 




rur^’ip.) . . ut^\p^) 

According to (6.6) pg. 332], the evaluation map satisfies the symmetry relation 


\p^) 


\p>^) 


and so we have 


LHS(|3) = ^ p,{y„ 2/„_i)V'^/,(g, Oil'll) . 

"*^0 {Pk-} K. 


(B.8) 


(B.9) 


The fundamental recurrence (© allows the summation over k in this expression to be performed, leaving 
us with 


LHS(^) = 


u^^\P^) 




(B.IO) 


But according to the definitions of and we have 

..., tj/n-l, 1)) = u^J^\P^). 

Substituting this in ( B.10|) then using the symmetry relation (B.8) in the form 

U^"\P^) ^ U^k\Pu) 

ut'\Pk) 4"\P^) 


gives the RHS of (B.7). 
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Appendix C 


Matrix Bessel functions 


Guhr and Kohler defined the matrix Bessel function 

= J dfi{U) exp(iTrt/1'x[/fc) 

where U G U{N;l3) with C7(iV;l) = 0{N) (orthogonal group), U{N;2) = U{N) (unitary group) and 
U{N;4) = Sp{N) (unitary symplectic group). Furthermore x = diag(a::i,..., xw), k = diag(A:i,...,/ cat) 
for P = 1,2 while x = diag(xi, xi,..., xat, xat), k = diag(A:i, fci...,/ cat,/cat) for /? = 4 (and thus each 
eigenvalue is doubly degenerate in this case). Let Un denote the Nth column of U (for /3 = 4 each element 
is itself a 2 X 2 matrix representing a real quaternion) and construct the random corank 1 projection of 
X by 

nxB, I1-=1 -UnUI^ (C.l) 

where 1 denotes the identity matrix. In the non-zero eigenvalues of this projection are termed the 
radial Gelfand-Tzeltin coordinates. 

Let A denote the N x N — 1 matrix with I’s down its diagonal and O’s elsewhere. Then we can write 
^]^\x,k) = J dfi{U) exp (^iTi xUkAA^^ exp (j,TrU^xUk{l — AA^)'^ . 

Noting that kA = Ak with k = diag(A:i,..., /cat-i) while fc(l — AA^) = (1 — AA^)kN, and defining x' as 
the {N — 1) X {N — 1) matrix 

x' := A^UAUA (G.2) 


we see that we have 


$ 


{x,k) = J dp{U) exp (^iTrix'k^ exp (^ik^ {Tr x — Ti x')^ . 


A simple calculation (see [g2|) shows that the eigenvalues of (C^) coincide with the non-zero eigenvalues 
of ( p.ll ). Gonsequently the second exponential in the integrand only depends on U^, so if we decompose 
dfi{U) as 

dfi{U) = dfj,{U)dUN 

where U G U{N — 1; /3) we have the factorization 


<Hx,fc) = 


dU^ exp (^ikj^{Tix — Trx')^ J dp{U) exp ^iTrix'fc^ 
AUn exp ^iA:A/(Trx — Trx^)^$|^lj^(x, k). 


Finally, the change of variables from Gat to the eigenvalues of (C.l), with the x’s given, is carried out 
according to Corollary |^. 


Appendix D 


Sampling from Jg"^’ ^ in a Monte Carlo calculation 

The eigenvalue PDF for Dyson’s circular ensembles of random unitary matrices with orthogonal (/3 = 1), 
unitary (/? = 2) and symplectic (/3 = 4) symmetry is given by 

n 

l<3<k<N 


|/3_ 


(D.l) 
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X 

Exact 

Monte Carlo 

0.2 

0.01687 

0.01690 

0.4 

0.2059 

0.207 

0.6 

0.6641 

0.674 

0.8 

1.1173 

1.15 

1.0 

1.2257 

1.30 

1.2 

1.0208 

1.15 

1.4 

0.8308 

1.05 


Table 1: The value of P 2 {x) in the case /3 = 4 as computed from (D.4) using M = 22, 500 distinct 
samples of 7 ^, tabulated against the exact value computed from ( p.5|) . The deterioration of the 
accuracy of the former as x increases is due in part to amplification of the error caused by the 
factor of in (D.4). 


The corresponding scaled two point correlation function is defined as 

~ 1) M _ g2ii-irc/V|/3 


P 2 {x) := lim 

N^oo Oat 


N 


X de^--- d6>Ar |1 - e*®' - e*®' 


iSi |/3|g2'n-ia:/V _ |/3 

1—3 3<j</c<A^ 


gi6lfc _ gi6lj 1/3 


(D.2) 


(note that with this scaling the mean eigenvalue spacing is unity). For general even /3 > 0 it was shown 
in |l^ that 

P2 {x) = {P/2f^^^0^{2Trxf(^cos27rx{J2xj - (^/^)) (D.3) 




where the average is over the PDF Jq^’~\x) ( ^.26 ) with k = 2//3, a = /3 = — 1 + 2//3. But we know the 
polynomial {x) computed from the random recurrence ([f.38|) with k = 2/P, ao = Po = —1 + 2/P has 
zeros with this distribution, and so X]j=i 2^3 = =: 7/3 where the symbol [a/*’] denotes the 

coefficient of x^. According to the general theory of Monte Carlo integration, if we sample 7^3 a total of 
M times we then have 




i=i 


1 


It is illustrative to compute (D.4) in the case /3 = 4, since then we have an alternative formula to 
(p.3|) for the exact evaluation 


, , , /sin27ra;\2 1 d /sin27ra; 


2tt dx\ X 


sint 


dt, 


(D.5) 


from which P 2 {x) can be computed to essentially arbitrary precision, thereby allowing the accuracy of 
(D.4) to be assertained. The results are listed in Table |^. 
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Appendix E 

Eigenvalue PDF for a generalization of the matrix structure ( p.5|) 

Here we will compute the eigenvalue PDF for the random matrix 


M = A + XBX^ + VtY 


(E.l) 


where A (B) is a fixed nxn (m x m) Hermitian matrix with distinct eigenvalue Oi,..., a„ (bi,..., b^), X 
is a n X 771 complex Gaussian matrix with real and imaginary parts having variance 1/2 and P is a random 
Hermitian matrix from the nxn GUE (see e.g. for the definition of the latter). Our expression will 
involve the operator 

^■'[/W](x) = lim [ e-+-^/V(s) xeR, 

jRe{s)=0 2771 

which is the inverse of the two sided Laplace transform 

/ OO 

e-^y{x)dx. 

-OO 


We will show that in the special case that t = 0 and B has rank 1 the term involving C ^ can be evaluated 
explicitly, reclaiming (^) in the case Si = 1 {i = 1,... ,n). 

Theorem 6. The eigenvalue PDF of the random matrix i[E.l\) is given by 


^ det 
n\ 


C det(l + Hs) ^]{xi — Oj) 


det[xl 




J det[a/ 


(E.2) 


Proof. Glearly, applying a unitary change of basis to A leaves the eigenvalue distribution of M un¬ 
changed; we may thus consider the random matrix 


M' := UAU^ + XBX^ + VtY 


for U Haar-distributed from U(n). Denote by E the operation of averaging over random variables, and 
consider the moment generating function (multivariate Laplace transform) 

fM'. E(exp(-Tr(M'C))) 

for Hermitian matrices C. We must therefore compute 

fM'{C) = E(exp(-Tr(t/H[/^C)))E(exp(-Tr(XBXl'C')))E(exp(-Tr(VtTG))) 

Now, if C has eigenvalues ci,..., c^, then 

m n m n - n 

E(exp(-Tr(NHXtC'))) = n n E(exp(-fo*G )) = n n = n d®t(l + He,)-' 

• 1-1 • 1 • 1 ■ 1 

and 

n n 

E(exp(-Tr((VtTG))) = E(exp(- y^((Vty,^ci)) = exp(-tc//2). 

Finally, by the Harish-Ghandra/Itzykson-Zuber formula (see e.g. |^ ) 

E(exp(-Tr(t7HDtC))) = (-l)"(-')/2 ^(i - 1)! det[exp(-a,c,)]y „ - 

det[a( J,j=i,...,ndet[c) 
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so we have 


2=1 


det[exp(—a^-Cj) exp(—ic^/2) det(l + Bci) 


(E.4) 


Note that this is analytic for C in a neighbourhood of 0. 

On the other hand, suppose we know that a random matrix M has distribution invariant under 
unitary change of basis, and has eigenvalue PDF 


1 


det[x^ det (xj)]i j'— 1 ^... 


(E.5) 


for some functions gj with Laplace transform defined in a neighbourhood of 0. Then 

E(exp(-Tr(CM))) 

= E(Eyec/(n)(exp(-Tr(CC/MC/t)))) 

=-—-— ■ - - dxi--- dxn det[exp{-XjC,)]ij^i^,,,^ridet[gj{xi)]ij=i^, 

Zdetlc^ | 2,7 = l,...,n J-OO J-oo 


Zdet[cr’ 


?-ll 


^det[c^ = 


Z det [c^ 


-ii 


dxi 


dxi 


dxn det 


' —OO 
/•OO 


e ^^*gj{x)dx 




dXn dot \£> ( 3 ^)] ^j=l ^ 


(E.6) 


Ji j'=l,...,ra 


where to obtain the second equality, use has been made of the Harish-Chandra/Itzykson-Zuber formula 
( |E.3| ). Comparing ( |E.6| ) with ( |E.3| ) allows us to deduce the value of Z and g[x), which when substituted 
in (E.5) gives (E.2). □ 


In the special case t = 0, B = diag(&,0,... ,0), the random matrix in (E.2) coincides with (p.5|). In 
this case 

£-ijgtsV2 Bs)~^]{xi - Qj) = £“^[(1 + bs)~^]{xi - aj) = 6 > 0 


and so 


det 


C det(l + Ss) ^]{xi — aj) 


J 2 ,j = l,...,n 

For oi > 02 > • • • > a„ the latter determinant is non-zero only when 

Xl > ai > X2 > ■ ■ ■ Xn > an 




in which case it is 1. Thus we see that (E.2) reproduces the eigenvalue PDF for the matrices (5.5) in the 
case that x is complex and the eigenvalues of A distinct. 

Another special case of interest is i? = 0. Noting that 

- aj) = (27rt)-i/2e-(=^^-“i)V2t 

we see from Theorem ^ that A -|- ^/tY for Y an element of the CUE has eigenvalue density 

1 , r det[x^ \i,j=l,...,n 






n!(27rt)"/2 -- -- det[af 

This is a well known consequence of the Harish-Chandra/Itzykson-Zuber formula (see e.g. p^). 
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